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Abstract We study the one-dimensional Burgers equation in the inviscid limit for Brownian initial 
velocity (i.e. the initial velocity is a two-sided Brownian motion that starts from the origin a; = 0). We 
obtain the one-point distribution of the velocity field in closed analytical form. In the limit where we 
are far from the origin, we also obtain the two-point and higher-order distributions. We show how they 
factorize and recover the statistical invariance through translations for the distributions of velocity 
increments and Lagrangian increments. We also derive the velocity structure functions and we recover 
the bifractality of the inverse Lagrangian map. Then, for the case where the initial density is uniform, 
we obtain the distribution of the density field and its n-point correlations. In the same limit, we derive 
the n— point distributions of the Lagrangian displacement field and the properties of shocks. We note 
that both the stable-clustering ansatz and the Press-Schechter mass function, that are widely used in 
the cosmological context, happen to be exact for this one-dimensional version of the adhesion model. 

Keywords Inviscid Burgers equation • Turbulence • Cosmology: large-scale structure of the universe 



1 Introduction 

The Burgers equation [lol | is a very popular nonlinear evolution equation that appears in many physical 
problems, see [6| for a recent review. It was first introduced as a simplified model of fluid turbulence, 
as it shares the same hydrodynamical (advective) nonlinearity and several conservation laws with the 
Navier-Stokes equation. Even though it was shown later on by [2§| and [13] that it can be explicitly 
integrated and lacks the chaotic character associated with actual turbulence, it still retains much 
interest for hydrodynamical studies. In particular, it can serve as a useful benchmark to test various 
approximation schemes devised for turbulence studies, since the nonlinearity is the same for both 
dynamics [l7l|. On the other hand, it has appeared in other physical situations, such as the propagation 
of nonlinear acoustic waves in non-dispersive media [23l |. the study of disordered systems and pinned 
manifolds [l6|, or the formation of large-scale structures in cosmology (26l ISlj . There, in the limit of 
vanishing viscosity, it is known as the "adhesion model" and it provides a good description of the large- 
scale filamentary structure of the cosmic web [s^l- In this context, one is interested in the statistical 
properties of the dynamics, starting with random Gaussian initial conditions [soL [27l | (i.e. "decaying 
Burgers turbulence" in the hydrodynamical context). Moreover, in addition to the velocity field, one 
is also interested in the properties of the density field generated by this dynamics, starting with an 
initial uniform density. 

This problem has led to many studies, focusing on power-law initial spectra (fractional Brownian 
motion) , especially for the twopeculiar cases of white- noise initial velocity [13, [13, El, [ll| or Brownian 
motion initial velocity [U, [H, [3| . The initial velocity fiuctuations are dominated by short wavelengths 
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in the former case and by large wavelengths in the latter case. In the present Universe, where the power 
spectrum is not a power law and converges at both ends, the velocity fluctuations are governed by 
scales that are somewhat larger than those where structures have already formed (thus the variance of 
the velocity field is still set by the linear theory) and this scale ratio was larger in the past (as the size 
of nonlinear structures was smaller) . In this sense the case of Brownian initial conditions is closer to the 
cosmological scenario. From the viewpoint of hydrodynamics, this is also an interesting configuration 
since in many hydrodynamical systems the power is generated by the larger scales. For instance, the 
Kolmogorov spectrum of turbulence, E{k) cx k~^^^, displays such an infrared divergence. Thus, the 
case of Brownian initial velocity was recently used in 20] to address the issue of local homogeneity 

In this article, we revisit the one-dimensional Burgers dynamics with two-sided Brownian initial 
velocity. In the spirit of the approach of [Tsj , using analysis methods (Laplace transforms) we obtain 
closed analytical results for n-point distributions (mostly in the limit where we are far from the ori- 
gin of the initial Brownian motion if n > 2). We check that our results agree with already known 
properties. In particular, we recover the property, derived by [§| through probabilistic tools for the 
one-sided Brownian initial velocity, that increments of the inverse Lagrangian map are independent 
and homogeneous. In our case this only holds for particles that are on the same side of the origin. 
We pay attention to issues that arise in the hydrodynamical context (e.g., velocity structure functions, 
Lagrangian displacement field) as well as the cosmological context (e.g., statistics of the density field, 
mass function of the collapsed structures associated with shocks). In particular, we compare our exact 
results with phenomenological models that are often used to describe the cosmological dynamics. 

We first describe in section [2] the initial Brownian conditions and the standard geometrical inter- 
pretation in terms of parabolas of the Hopf-Cole solution of the dynamics • Adapting to our case 
the method presented in [l^, this will allow us to express all statistical properties in terms of the 
transition kernel associated with Brownian particles moving above parabolic absorbing barriers. We 
present this propagator in sect. [H decomposed over a continuous set of eigenfunctions built from the 
Airy function (whereas the white-noise case leads to a discrete spectrum, that is also built from Airy 
functions). Then, we derive closed analytical expressions for the one-point velocity distribution Px(v) 
in sect. 21 as well as the distribution, Px{q), of the initial Lagrangian position q of the particle that 
is located at the position x at time t. Next, we study the two-point and higher-order distributions in 
sect. El and we obtain simple analytical results in the limit where all particles are far from the origin. 
This allows us to derive the distribution of the density field in sect. |6l for the case of a uniform initial 
density. Next, we consider the statistics of the Lagrangian displacement field in sect. [71 In the same 
limit where the particles are far from the origin, we obtain the n— point distributions, Pq^ixi), of the 
positions Xi at time t of the particles that were initially at positions g^. We also derive the probability 
-shock ^T^g^^ ^y^Q particles initially separated by a distance q have coalesced into a single shock by time 
t. Finally, we obtain in sect. |8l the mass function of shocks and their spatial distribution. 

The reader who is not interested in the technical details of our derivations may directly go to 
section O to survey most of our practical results. 



2 Initial conditions and geometrical solution 



We consider the one-dimensional Burgers equation for the velocity field v(x,t) in the limit of zero 
viscosity, 

dv dv d^v .., „j_ 



dt dx dx"^ 



with 



V 0+. 



As is well-known [29], [13|, introducing the velocity potential ^p{x,t) and making the change of variable 
'ijj{x,t) = —2vln9(x,t) transforms the nonlinear Burgers equation into the linear heat equation. This 
gives the explicit solution 

{x-qf Vo(g) 



v(x,t) — —- with ip(x,t) 
ox 



-2v\tl 



_dq_ 



exp 



2i/ 



(2) 



where we introduced the initial condition ipoil) — V^(9j^ — 0). Then, in the limit v ^ 0^ the steepest 
descent method gives 

{x - qf 



ip{x, t) = min 

9 



2t 



and v{x,t) = 



X - q{x, t) 



t 



(3) 
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where we introduced the Lagrangian coordinate q{x,t) defined by 

(x- q)^ . . . 

4'oi.q) H 2^ is minimum at the point q = q{x,t). (4) 

The Eulerian locations x where there are two sohitions q_ < q^ to the minimization problem ^ 
correspond to shocks (and all the matter initially between and (7+ is gathered at x) . The application 
q 1-^ x{q,t) is usually called the Lagrangian map, and x 1— > q(x,t) the inverse Lagrangian map (which 
is discontinuous at shock locations). For the case of Brownian initial velocity that we consider in this 
paper, it is known that the set of regular Lagrangia.n points has a Hausdorff dimension of 1/2 [43l|. 
whereas shock locations are dense in Eulerian space [ij, |4lj . 

In this article, we take for the initial velocity field vq (q) a bilateral Brownian motion starting from 
the origin vo{0) = 0, and we also normalize the potential ipo by "00(0) = 0. Thus, introducing a Gaussian 
white noise ^(g), we can express the initial conditions by 

vo{q)= r dq'aq'), ^0(9)= Tdg' r dg"e(g")- (5) 
Jo Jo Jo 

All initial fields are Gaussian and fully determined by their two-point correlation, which we normalize 

by 

(C(g)) = o, maq'))=D5{q-q'), (6) 

where (..) is the average over all realizations of ^. This gives for instance 



{voiqi)vo{q2)) = Dqi, (■0o('7i)V'o(g2)) = y 
and for the initial velocity distribution at location q, 



' 2 9i 

9192 - -IT 



for < gi < (72, (7) 



t = 0: p^{v) = ^^e-^'/^'<^ with cTliq) = Dq. (8) 



Note that the initial fields over the two sides g < and q > are independent. The initial velocity 
^0(9) is not homogeneous, since the origin q = clearly plays a special role, but it has homogeneous 
increments, as seen from the equality, 

ior any qi,q2 : ^0(92) - "0(91) = / dq£,{q), {[vo{q2) - vo{qi)f) = D\q2 ~ qi\. (9) 

J qi 

Then, the energy spectrum Eo{k) of the initial velocity field is 

Eo{k) = —k-\ with (K (92)- (91)]') =2 / dfc(l-e"=te'''^))i?o(fc). (10) 

Thanks to the scale invariance of the Brownian motion, the scaled initial potential ijjo{Xq) has the 
same probability distribution as X^^'^ipoiq), for any A > 0. Then, using the explicit solution ^ we 
obtain the scaling laws 

^{x,t)^= t^iP{x/t^,l), v{x,t)^= tv{x/t^,l), q{x,t)^l^ t\{x/t^,l), (11) 

where ^= means that both sides have the same probability distribution. Thus, any equal-time statistics 
at a given time t > can be expressed in terms of the same quantity at the time t = I through 
appropriate rescalings. In this article we only investigate equal-time statistics, so that t can be seen as 
a mere parameter in the explicit solution ^ from which we derive our results. 

In the cosmological context, the time t in the Burgers equation ([1]) actually stands for the linear 
growing mode D^{t) of the density fiuctuations, the spatial coordinate a: is a comoving coordinate 
(that follows the uniform Hubble expansion) and, up to a time-dependent factor, the velocity v is the 
peculiar velocity (where the Hubble expansion has been subtracted), see (2^, [HI. In these coordinates, 
the evolution of the density field is still given by the continuity equation (|106[) below, where the density 
p is the comoving density. If we take i/ = 0, that is we remove the right hand side in Eq.([T|), this is 



4 



the well-known Zeldovich approximation [52, l47l | , where particles always keep their initial velocity and 
merely follow straight trajectories. The diffusive term of ([T|) is then added as a phenomenological device 
to prevent particles from escaping to infinity after crossing each other and to mimic the gravitational 
trapping of particles within the potential wells formed by the overdensities [1^. Of course, this cannot 
describe the inner structure of collapsed objects (e.g., galaxies) but it provides a good description of 
the large-scale structure of the cosmic web [3J] . 

As is well-known [lol |. the minimization problem (j4|) has a nice geometrical solution. Indeed, let 
us consider the downwarcQ parabola Vx.cio) centered at x and of maximum c, i.e. of vertex {x, c), of 
equation 

^.,c(g) = -^^^ + c. (12) 

Then, starting from below with a large negative value of c, such that the parabola is everywhere well 

below ipo{q) (this is possible thanks to the scaling tpoiXq) ^= A'^/^'(/'o('7) which shows that ipoiq) only 
grows as Iq]^'"^ at large \q\), we increase c until the two curves touch one another. Then, the abscissa 
of the point of contact is the Lagrangian coordinate q{x,t) and the potential is given by 'il}{x,t) = c. 
(We show below in Fig. [T] the case where the Lagrangian coordinate q'{x,t) is somewhere in the range 
< g' < q.) 



3 Transition kernel with parabolic absorbing barrier 

For the Brownian initial conditions (O, the process q i— > {"00, ''^o} is Markovian, going from q ^ Q 
towards positive or negative values. Then, following the approach of 18] (where it was applied to 
white-noise initial velocity), from the geometrical construction (|12[) we can see that a key quantity 
is the conditional probability density ifa;,c(9i, V'l, ^^i; 92, V'2, ^2) for the Markov process {"00(9), i'o('7)}, 
starting from {ij^i^vi} at qi > 0, to end at {ip2,V2} at q2 > qi > 0, while staying above the parabolic 
barrier, ■00(9) > 'Px,c{q), for qi < q < 92- It obeys the advective-diffusion equation 

D 92 

K.j;^c{qi,ipi,vi;q2,'il^2,V2) = —^^Kx,c{qi,ipi,vi]q2,il)2,V2) (13) 
over the domain ^ > Vx.dq), with the initial condition at q2 — 91 

-?^x,c(9l, V'l, «i; 91, 02, W2) = ^(02 - '0l)^(w2 - Vi), (14) 

and the boundary condition 

Kx,c{qi,ipi,vi;q2,ip2,V2) =0 at 0^2 = ■Pi.c(92) for W2 > — 7^(92)- (15) 

dq 

Equation p3| is also the Klein-Kramers equation for the distribution function P{x,v;t) of Brownian 
particles, in the limit of zero external force and zero friction coefficient but finite diffusion coefficient, 
where we identify the position, velocity and time coordinates as {x, v; t} = {■02, V2', 92}- The boundary 
condition (fT5|) simply means that particles cannot come back from the absorbing region (i.e. curves 
that cross the parabola are "lost" and do not contribute to the probability density Kx^c)- 

In the case of white- noise initial velocity studied in ^Ts'l , the velocity potential ip itself is a Brownian 
motion so that the relevant propagator only involves one dependent variable, 0, as K^]^ {qi^ ipi] 52, "02)- 
In our case, since is now the integral of the Brownian motion w, the propagator Kx^c introduced 
in ()13|) involves the two dependent variables v and 0. Thus, we have a diffusion in a two-dimensional 
{0, w}— space rather than the one-dimensional 0^— space as in [Tsj . As we shall see below, the propagator 
Kx^c involves an expansion over a continuous spectrum of eigenfunctions that are built from the Airy 
function, whereas the white-noise case leads to a different expansion over eigenfunctions that are still 
built from the Airy function but form a discrete spectrum, see (isj. 

^ In the literature one usually defines the velocity potential as « = —dxtp, which leads to upward parabolas. 
Here we prefer to define v = dx^ to simplify the interpretation of the process (q, i/jq, «o) in terms of the dynamics 
of a Brownian particle. 



_d_ _d_ 

dq2 902 
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The conditional probability density Kx^c associated with the left-handed Brownian motion q2 < 
gi < can be obtained from the symmetry q — > —q as: 



< 91 < 92 : Kx,ci-qi,'4'i,vi;-q2,ip2,V2) = ^'"i'^ 12,'4'2, -V2), 



(16) 



hence we only need consider Eq. (|13p for < gi < 92- To solve this equation it is convenient to make 
the change of variables 



KxA<li^'4'i^Vi]q2,i>2,V2) = lC{qi,yi,Wi]q2,y2,W2), 

dVx,c 



with y ~ VxAl) = H ^ c, w 



dq 



{q) = v + 



q — X 



(17) 
(18) 



to obtain a simpler boundary at the fixed vertical half-line {y = 0,w > 0) in the (y, w) half-plane for 
/C, J/ > and —00 < w < 00, instead of the parabolic boundary for K. From Eq. (|13p the kernel /C 
satisfies the equation with constant external force 



_d_ 

dq2 



d 



1 d 



W2 



dy2 t dw2 



Then, making the transformation 



/C((ji, 2/1, wi; (72, 2/2,^2) = — G(r;ri,Mi;r2,M2) exp 



W2 -wi 92 - 91 



Dt 



2Dt^ 



with 



T = q2-qi, r=\l — y, u=^ — w, 



we obtain the simpler advective-diffusion equation for r > and a; > 0, 

dG_ _ d^G 

dr ^'^ dr2 du\ ' 

with the initial and boundary conditions 

G(0;ri,ui;r2,U2) = ^(»'2 - ?'i)'5('U2 - ui), G(t; ri, mi; 0, M2) = for U2 > 0. 



(19) 

(20) 
(21) 

(22) 
(23) 



Thus, G(r; ri, mi; r2, 1*2) is the conditional probability density of Browni an p articles with unit diffusion 
coefficient and absorbing barrier at r = 0. This quantity was obtained in [11| and we briefly recall below 
his procedure using our notations. We first take the Laplace transform of G as 



G(s;ri,ui;r2,M2) 



dre ''^G(r;ri,Mi;r2,U2), 



hence Eq. (|22)l gives 



d d"^ 
dr2 du2 



G(s;ri,Mi;r2,U2) = 5{r2 - ri)5{u2 - ui). 



(24) 



(25) 



Next, to obtain an ordinary differential equation, it is convenient to expand over the eigenfunctions 
^-f ^^gg ,^{u2) associated with Schrodinger's equation 



dii2; 



gsM{u) = 0, whence Qs^iu) = Ai 



-i^u H — - 



(26) 



using the fact that the standard Airy function Ai(a;) is the only solution of Ai"(x) = xAi(x) that 
vanishes at both ends x ±00 [ij]. We recall in Appendix lAl some useful properties of this entire 
function. Using the integral representation (|165p . we obtain the orthogonality property 



/•oo 




s 


Ai 




S ' 


/ du u Ai 






— V u - 




y —oo 













(27) 
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and the closure relation 







S ' 


Ai 




S ' 


J —oo 


—vu - 




—vu' - 





■S{u-u'). 



Therefore, we can see from Eas. (P5)) - (P51) that Eq. ips)) has the particular solution 



G'o(s;ri,ui;r2,U2) 



(28) 



-'■i^Si^Ai 


—vui ^ 


S ' 


Ai 


s ' 










L v"^ J 


- ra) + 0{ 




-ri) 







(29) 



where 6 is the Heaviside function. We can check that Go vanishes for \r\ oo and for |u| oo. Then, 
since we have not taken into account the boundary condition at r2 = of (|23p yet, Go is the Laplace 
transform of the probability density of Brownian particles over the unbounded plane (r, u) (thus Go 
only depends on the length \r2 —ri\). Note that the solution to this unbounded problem is well known 
to be the Gaussian JJj] 



Go(T;ri,ui;r2,U2) 



V3_ 



r('"2-''l-MlT)^ + -%(r2-ri-«iT)(ll2-til)-i(ti2-Ml)^ 



(30) 



as can be checked by substitution into Eg. ([2211 . Therefore, Eq. ([30|) is the inverse Laplace transform of 
Eq.dMl). 

Next, in order to satisfy the second constraint ()23p . we must subtract to Go an appropriate solution 
Gi of the homogeneous form of Eg. ([25]) . From Eg. ([26|) . we can see that Gi can be written as a 
combination of eigenfunctions e^'' gs,p,iu2), that must be restricted to /i > to ensure that G 
vanishes for r2 +oo. Moreover, for r2 = only the first part 9{—v)9{ri — r2) contributes to Go in 
Eq. ((29| . Therefore, to compensate for this term at r2 = for U2 > 0, we must look for a function Gi 
of the form 



Gi(s;ri,ui;r2,U2) 



dh' e 



3i/Ai 



(t)sM{r2,U2), 



where the function (j)s^^(r,u) can be written as 



Ai 



-nu H TT 



with some weight Ws^^{fJ.), and satisfies the constraint 



(t>s,,y{r = 0, u) = Ai 



for u > 0. 



(31) 



(32) 



(33) 



This is a half-range problem as we must decompose a given function (here Ailvu + s/v'^]) over half the 
domain [u > 0) using only half of the eigenfunctions gs,fj,{u). Using the results of j31| . who studied the 
Klein-Kramers eguation, and taking the limit of zero friction but non-zero diffusion, [Tlj obtained: 



1/ > : 4>s,u{r,u) 



d/i 3l/-^/^U^/^ _2„3/2/ -3 , -3> _,3 



-^U -I 77 



(34) 



27r v"^ + fj,'^ 

Substituting into Eq. pTj) . we obtain for the solution G of Eq. (|^5|) . with the boundary conditions (pS)) . 

G = Go-Gi, with (35) 



Gi = 



di^d^ 9i/^/^^^/^ _2„3/2^,^-3j_„-3,| _„3„ 3, 

27r v^ + 



■ e 3' 



S ' 


Ai 


s 


vul H 


-M"2 H 7 


. 







(36) 



We describe in Appendix [B] how the solution can be directly obtained for the half-range expansion 
problem ([32|) - (j33|) . associated with the Brownian dynamics ([22l) . rather than first solving the problem 
associated with the Klein-Kramers dynamics and next taking the limit of zero friction, see Eg. (|188p . 
This also allows us to derive the more general identities (|187p . (|189p . that we need in the following 
sections. 
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We can see from the explicit expressions (|29|) . (p6| . that the kernel G also satisfies the backward 
evolution equations (compare with Eqs. ([22ll . ([25| ) 



d 




ori 


dul 


d 




Ul^T— 


dul 



G(r;ri,ui;r2,U2) = 0, 

G{s\ ri, ui; r2, U2) = ^(^2 - 7-i)(5(m2 - wi), 



as well as the boundary condition (compare with (I23p ) 

G(t; 0, ui; r2, U2) = for wi < 0. 



(37) 
(38) 

(39) 



Equation (f39|) merely states that the trajectory r(T) starting on the absorbing barrier at ri = must 
start in the upward direction mi > not to be immediately absorbed. 

For later calculations we also need two kernels A and H that are derived from G. Thus, we define 
the propagator A, that will be associated with Brownian particles that come within a small distance 
e from the parabolic absorbing barrier, by 



Z\(T;ri,Mi;r2,-U2) = hm -[G{T]ri +e,Mi;r2 + £,1*2) - G'(r; ri, wi; r2, U2) 



(40) 



From Eas. ((29|) and ([36]) we have for its Laplace transform A 



Zi(s;ri,ui;r2,M2) 



27r 



g^3/2^3/2 g-§s-'/^(^-^+p-3) g-^^ri-M-^rs 



xAi 



s 


Ai 


s 


Z^Ml H — ^ 


~^U2 + 


L i/^ J 







(41) 



Next, we define the kernel Hoc{ri,ui), associated with Brownian particles that stay forever above the 
parabolic absorbing barrier, by 



i?oo(?^i, ui) = lim e '^''^ H{T;ri,ui), 

T — ^ + C30 
/•oo nOG 

with H{T;ri,ui)~ / dr2 / du2 e^^^''' G'(t; ri, ui; r2, M2). 

Jo J -00 



(42) 
(43) 



Using Eas. ([29|) and ([36]), and the property (|172p . we obtain after integration over ^2 and U2 for the 
Laplace transform H, 



H{s;ri,ui) 



di/— Ai 



-VUl + 



3^" 



27r + 



iy)-6'(-i^)(l-e'^ '■I) 
Ai 



t^Ul H ^ 



g(#-5^)/f'^ (44) 



The behavior for t ^ cx) of H{T;ri,ui) is determined by the rightmost singularity of H, which is 
located at s = I/72. At this point, the first integral in Eq. ([44)) diverges for 1/ ^ 0+ whereas the second 
integral diverges for /i — 0+. Therefore, the singularity is governed by the behavior of the integrand 
for 1/^0+ and /i 0+, so that we can expand the first Airy function and the ratio l/(i^'^ + /i'^), which 
yields 



7 



H 



1 



,ui/-f 



s — 7 

This gives for the function _ffoo(''i7 wi): 



3^-3/2g-f.--.-Vi/7^Ai 



U^ 1 

7 1/"^ 



Hoo{ri,ui) = e 



— p"i/7 



^3:.-3/2e-i''-^-''^-/^^Ai 



7 



(45) 



(46) 
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Finally, using the transformations (fT7|) and (f20|) , we obtain in terms of the original variables 

if,,,((7i, 1^1, t;i; g2>2, ^^2) d^-sdi-a = h'" ^'^^^-^^) h n, mi; r2, ua) drsdus, (47) 

with 



T = 7'(Q2-Qi), = 27^ 



l2 



, = 27(K, + Q,-X). (48) 



Here we introduced the dimensionless spatial coordinates (which we shall note by capital letters in this 
article) 



y = — = 777^7, whence X = Q + y for regular points. (50) 



and the dimensionless velocity 

tv V 

In a similar fashion, the dimensionless velocity potential coordinates in (|48p are 

f'^^ and C=^. (51) 

Next, from Eq. (l42p the kernel associated with Brownian particles that remain forever above the 
parabola Vx,c reads as 

lim / dV'2dw2i^x,c('7i,'0i,'yi;92,^2,W2) = e""i/''i7oo(?'i,Ui), (52) 

92— '+00 J 

whereas the propagator associated with Brownian particles that come within a small distance dc from 
the parabolic absorbing barrier is from Eg. ((40)) 

lim -^[Kx^ciqi,ipi,vi;q2,ip2,V2) ~ i^a;,c+5c(gi, V'l, ^^i; 92, "02, ^2)] dV'2du2 = 

(5c— fO be 

1- e'-h^+{n2-u,)/^ r^^ui-r2,U2) Ar2Au2. (53) 
7 

4 One-point distributions 

4.1 Results for arbitrary Eulerian location x 

In this section we consider the one-point velocity distribution Px{v) at the Eulerian location x. From 
the explicit solution (jS]) , it can be derived from the probability distribution px (<?) of the Lagrangian 
coordinate q{x,t). Thus, we have from Eq.(l3]) 

Px{v)^tpx{q) and q = x - vt, (54) 

where we note px (v) and Px {q) the probability distributions of the velocity v and of the Lagrangian 
coordinate q, at the Eulerian location x and time t. Here we used the property that q{x, t ) is well 
defined for any x except over a set of zero measure in Eulerian space associated with shocks [4l| . 

Then, from the geometrical construction (fT2|) . we are led to consider the bivariate probability 
distribution, px(0 < 9' < 9; c)dc, that the first contact point of the potential ipoiq') with the family 
of downward parabolas Vx,c{q')i with c increasing from —00, occurs at an abscissa q' in the range 
< g' < (7, with a parabola of height between c and c + dc. This will give us in turn the cumulative 
distribution Px(0 l£ q' l£ q) by integrating over c. Then, for g > 0, we can write this probability 
distribution as 



Px{0 ^ q' < q,c)dc — lim / d7/'_di>_d^di;d^+du+ i^Tj; c(0, 0, 0; V"-, i'-) 

q±^±oo J 

X [Kx,c{0, 0, 0; q, i/j, v) - Kx.c+ddO, 0, 0; q, ip, v)] Kx,c{q, 4', u; g+, 1/'+, i)+), (55) 
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Fig. 1 (color online) Geometrical interpretation of the initial conditions tpo{q") associated with the probability 
Px{0 < q' < q, c)dc. The Brownian curve tpo{q") is everywhere above the parabola Vx.c and goes below 'Px,c+dc 
somewhere in the range < g" < q. From the constraints ipo{0) ~ and i/'o(0) = 0, see Eqs.([5}, it goes through 
the origin with an horizontal tangent. To obtain the cumulative probability, PxiO < '?'<'?), we must then 
integrate over the height c of the parabola. 



where we used the Markovian character of the process q i— *■ v}. Thus, we could factorize in Eq. (|55p 
the probability (0 < q' < q, c)dc into three terms, which correspond to the probabilities that i) tpoW) 
stays above Vx.c for q' < 0, ii) i^aiq') stays above 'Px,c, but does not everywhere remain above 'J^x,c-\-dc^ 
over the range < g' < g, while reaching an arbitrary value v} at q, over which we will integrate, 
and iii) V'o('z') stays above 'Px,c for q' > q- We show in Fig. [l]the geometrical interpretation of Eq. ([55|) 
(where we did not try to draw an actual Brownian curve ^poil) which has no finite second-derivative). 

We can easily check that in the limit x — *■ -|-c» and q — *■ +00 with q ^ x, the integral over c of 
Eq. lfSSj) gives unity as it should. It is convenient to first compute cumulative probabilities as in ([55]) and 
to take the derivatives afterwards to derive the probability densities. This ensures that probabilities 
are well normalized and it avoids coming across ill-defined expressions. Indeed, since the curve ipoiq) 
has a continuous derivative, it is tangent to the parabola Vx,c at the first contact point. Then, this 
point corresponds to r = and u = in terms of the reduced variables (j2ip . where the Brownian 
kernels are singular. For instance, the expression is not well defined if we naively put ri = r2 = 0. 
Other ambiguities or seemingly divergent quantities appear if we try to directly compute probability 
densities by using Taylor expansions. 

Then, using the relations (|T7)) - (|53p . we obtain 

Px{0 < q < q,rQ)dro = e^'^/^ dro j drduHao{rQ,u)A{q;r^,-u;r,u)H^{r,u), (56) 
where we defined 

Using the results of section [3l the integration over r and ro gives 

:^e(^-i)Q/(s)J(s,2X), (58) 
-ioo 27ri 

where we introduced the dimensionless variables Q and X as in and we defined the functions 



/C30 
dz J(s,2z), 
-00 



(59) 
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For y > 0, we obtain using Eqs. p87|) - p90| . 
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For y < 0, using Eq. (|189p in the second term of Eg. ([SO)) , we obtain 



y<0: J(s,y) = 6 / 
Jo 



dA* . 3/2^^. 



Ai 



1 ■ 


Ai 


s 


-m + — 


-My + — 









(60) 



(61) 



(62) 



Therefore, since we have the primitive 



Au Ai 



Ai 



S2 



{A. 



Sl - S2 





Sl" 


Ai 


i/U - 


S2' 


- Ai 


z^w - 


Sl" 


Ai' 


i^u - 


S2" 






















" J/2. 



}, (63) 



the integral (15911 reads as 

3 



di/ 



g§(.^/=-i).- 



-Ai 



Using Eq. (|199p this yields the simple result 



/(.) = 



s - 1 



Therefore, in terms of dimensionless variables, the cumulative probability (|58|) reads as 

J^too 27ri s - 1 



(64) 



(65) 



(66) 



On the other hand, since the system is statistically invariant through reflection about the origin, we 
have the symmetry pa;(0 < q' < q) ~ p-x{^q < < 0). This implies that the cumulative probability 
distribution associated with a Lagrangian coordinate q' on the negative real axis reads as 



+'°° ds 



iiQ J{s,-'2.X) 



s- 1 



Px(-g<g'<o)= / -^e(^"i)'? 

We can also check Eg . ([57)) through an explicit calculation similar to ([55)1 . 

From Eq. ([M)) and Eqs. ([55|) - ([57)l . the cumulative velocity distribution is given by 

w < ^ : Px{v <v' < x/t) =P:r(0 <q' <x-vt)^ Px{0 <Q' <X -V), 

w > I : Px{x/t <v' <v) ^px{x -vt<q'<0) = P_x(0 < Q' <V - X), 

where we introduced the dimensionless velocity V defined as in ([50)1 . Of course, Egs. ([00)l - ([00)l agree 
with the scalings ()lip . 

Letting \Q\ — > oo in Egs. ([66)) - ([67| . or \V\ — > oo in Egs. ([68)) - ([69)l . we obtain the probabilities that 
the Lagrangian coordinate q, associated with the Eulerian coordinate x, is located on either side of the 
origin (or that the velocity v is smaller or greater than x/t): 



(67) 



(68) 
(69) 



Pxiq > 0) ^ px{v < x/t) = J(l, 2X), 
Px{q < 0) = px{v > x/t) = J(l, -2X). 



(70) 
(71) 
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V=v/(aDt) log[|V[]=log[|v|/(aDt)] 



Fig. 2 (color online) Left pan el: The probability distribution Po{V) of the reduce d v elocit y V — v/{2Dt) at 
the origin x — 0, from Eq.((74}. The dashed lines show the asymptotic behaviors l|76p and (|77p . Right panel: 
Same as left panel but on a logarithmic scale. 



Here we used the fact that the large-Q behavior of Eqs. ([66|l -([67 l) is set by the rightmost singularity of 
the ratio J(s, ±2X)/{s— 1), which is the simple pole at s = 1. From Eg s. ([61 ]) - ([62)1 we obtain for x > 0: 



x>0: Pxiq<0) 



di/ — Ai 



v2X + — 



Px('7>0) = l-p,(q<0). 



(72) 



We can check that the sum of these two probabilities is equal to unity. As expected, Eg. ([7^ shows 
that pxisi < 0) decreases as x gets larger and it goes to zero for x — > +oo. For a; = 0, both quantities 
are equal to J(l, 0) = 1/2, as can be checked from the explicit computation of the integral in Eq. (TfS]) . 
Finally, from Egs. ((66| - ([67l) the probability densities are given by: 

^e^^-''>Qj{s,2X), Px{-Q)^ ^e(^-i)«J(s,-2X). (73) 

-ioo 27ri 27ri 



This also gives the velocity distributions through the relation ((50 



4.2 Velocity distribution at the origin a; = 

We consider here the one-point distribution at the origin a; = 0. From Eg. ([73]) we can check that the 
distribution is even (we have Q = ~V at X — 0). For y > it is given by 

V>0: Po(V^)= r'":^e(^-i)^J(s,0), PohV) ^ Po{V). (74) 
J-joo 27ri 

The behavior for 0"*" is determined by the behavior at s ^ +oo of J(s,0). From Eg. ([62)l and 

using Eg. (|177p we obtain 

s^+oo: J(s,0) - — (75) 

TT 

which leads to 

V~.0+: PoiV)^-^^^. (76) 

TT V TTV 

Thus, we obtain an inverse square-root divergence for Po{V) at 0. 

The behavior of Po{V) for large V is governed by the rightmost singularity of J(s,0), located 
at s = (associated with the branch cut along the negative real axis). There, J(s,0) behaves as 
J(s, 0) ~ s~^/^, because of the first term in Eg. lpTj) . This yields 

V-^+oo: Po(V^)^ j^^^-'/'e-^. (77) 
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Note that initially, at time t = 0, the velocity at the origin is not random as it is equal to zero, see 
(l5|). Then, for t > the nonlinear evolution of the velocity field v{x, t) broadens this initial Dirac peak 
and gives rise to the exponential tail (ffT]) at large velocities and to the power-law peak (|76p at low 
velocities. We show in Fig. [2] the velocity distribution Po{V), as well as the asymptotic behaviors (|76p 
and (|77p. that happen to describe very well most of the distribution. 

We can note that since all quantities can be expressed in terms of the scaling variables p9|) - ([50)) . the 
exponential tail (|77p can be understood from simple scaling arguments applied to the initial velocity 
field. Thus, for a particle of initial Lagrangian position g > to reach the Eulerian position x = at 
time t, we can expect its initial velocity to be of order vq ^ —q/t. From Eq.® this corresponds to a 
probability of order e^^°^^^''^o^'^'>^ ~ e""^/*^ ~ e^'^ , where we did not write factors of order unity in the 
exponent, which cannot be obtained by such arguments. Thus, we recover the exponential tail (|77p (at 
X = we have V = -Q). 



4.3 Velocity distribution for |a;| oo 

Finally, we consider the one-point velocity distribution at large \x\. By symmetry, we only need consider 
X +00. Using the relation X — Q + V and Eg. ([75)1 . we can write the velocity distribution in terms 
of the reduced variables X and V as 

— e(^-i)(^-^) J(s,2X). (78) 
-ioo 27ri 

We now consider the limit X +00 at fixed V. Then, making the change of variable s = 1 + ifc, we 
obtain at leading order (fc being of order 



00 



Px{V) ^ I ^giHx-v) ^{i-VT+ik)2X ^ ^ g-ifcy-xfeV4 ^ f ^ (79^ 



2vr 7_oo 27r 

Therefore, in terms of the variable u, we recover as expected the initial Gaussian ([8]). This can be 
understood as follows. A remote region [x — L/2,x + L/2], with L ^ x, has a mean initial velocity 
Wo ~ V Dx that is much larger than its initial velocity dispersion Avq ^ \J DL, see Eq.®. Then, 
this domain remains well-defined and not strongly disturbed by neighboring regions until times of 
order with Avoti, = L, that is Dt^ = L. Conversely, at any time t, for x :s> Dt^ (i.e. X ^ 1) 
it is possible to make such a separation of scales and to identify a region of size L around x, with 
Dt^ L X, that moves in a collective fashion with a mean velocity ~ vq{x) that is set by the initial 
velocit y. T herefore, we recov er at leading order the initial Gaussian velocity distribution, of variance 
(Ty = V Dx (i.e. cry = \/X/2 3> 1), and the nonlinear evolution only modifies the velocity distribution 
by changes of order Av ^ Dt (i.e. AV ^ 1). The result ([79]) confirms this simple scaling argument. 
This is an illustration of the "principle of permanence of large eddies" [27,] , that holds for more general 
energy spectra, Ep(k) oc fc", with n < 1. As suggested by this discussion, and as checked in numerical 
simulations [3, US] , the stability of large-scale structures is not only a statistical property but actually 
holds on an individual basis, that is for each random realization of the velocity field. 

Of course, this reasoning does not apply to rare events, such as those where the displacement x — q 
remains of order x. In particular, from Eq. (|72p . we obtain for the cumulative probability to have a 
negative Lagrangian coordinate q the asymptotic behavior 



-|-oo : 



P'Aq <0)^Px{v> x/t) - {s^T^/iX^ e-^^^ . (80) 



Thus, we obtain an exponential tail for these very rare events. It can again be understood from simple 
scaling arguments, as for the exponential tail ([77]) . Thus, for a particle with Lagrangian coordinate 
(? < to reach the position x ^ 2Dt^, we can associate the initial velocity vq = {x — q)/t and the 
probability e~^^~''' using Eq.Q without writing numerical factors. Then, the maximum over 

g < of this exponential weight is reached for q = —x, which gives a weig ht ~ e"^/* ~ e"^ that 
agrees with Eg. ([50]) . We show in Fig. [3] the probability (9 < 0), as well as the asymptotic decay (|5n|) . 
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X = x/(2Dt=) log[X] = log[x/(aDt=)] 

Fig. 3 (color online) Left panel: The probability Px{<l < 0) (equal to the reduced cumulative probability 
Px{Q < 0)), that a part icle, located at a position a; > at time t, was i nitially located on the negative real 
axis, g < 0, from Ea. (|72p . The dashed line is the asymptotic behavior (|80ll . Right panel: Same as left panel but 
on a logarithmic scale. 
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Fig. 4 (color online) Geometrical interpretation of the initial conditions 4^o{q) associated with the contribution 
^^1,12(0 ^ 91 ^ Qi,ci;q2 > q2)dci. The Brownian curve tpo{q) is everywhere above the parabola Vxi^ci, it goes 
below 'Pxi,ci+dci somewhere in the range < gj < qi, and it goes below the parabola 1^x2,0,, of center X2, that 
intersects 'Pxi,ci at q, =52- This counts all paths with a first-contact parabola 1^x2,02 such that C2 < c, and 
g« > 52 (which implies q'2 > 52)- 

5 Two-point and higher-order distributions 

5.1 General results for xi < X2 and < gi < 52 

We now study the two-point Eulerian velocity distribution Pxi.x2 (^i: ^^2), with xi < X2- As in sectional 
we first consider the distribution Pxi,x2i'liT Q2) of the Lagrangian coordinates qi,q2, associated with 
the Eulerian positions xi,X2- For the Brownian initial conditions ©-([ni) shocks are dense [H, |4l|. 
Therefore, for xi < X2 there is almost surely a shock between xi and X2 and these two Eulerian points 
are associated with two different Lagrangian coordinates qi 7^ 92- This can also be understood from 
the fact that at the contact point qi (resp. (72) the curve ipoiq) is tangent to a parabola Vxi.ciiq) 
(resp.7^x2,c2 from the geometric construction recalled in (fT2|) . Then, since two parabolas 1^x1,01 
and 'Px2,c2 with xi 7^ X2 have different tangents at any point q (indeed d7-'x,c/d(j' = —{q — x) jt), the 
curve Ipoiq) cannot be tangent to both parabolas at a common point qi = q2 (in both steps we used 
the property that the derivative tpoiq) is continuous, being a Brownian motion). Therefore, we almost 
surely have qi ^ q2- Then, since particles do not cross each other we have qi < q2 for xi < X2- 

As in sccti0n r4.ll we first consider the cumulative probability distribution, Pxi,x2iO ^ 9i ^ <?! ! 92 ^ 
(72), that the Lagrangian coordinates q'i,q2, associated with the Eulerian positions a:i,a;2, are within 
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Fig. 5 (color online) Geometrical interpretation of the initial conditions ipo{<l) associated with the contribution 
0:2(0 ^ Qi ^ '?iiCi;g2 > 92, C2)dcidc2. The Brownian curve ipo{<l) is everywhere above the parabolas Vxi.ci 
and 'Px2,c2, it goes below 'Pxi,ci+dci somewhere in the range < q'l < qi, and below the parabola 'Px2,c2+dc2 
somewhere in the semi-infinite range q'2 > 92. The height C2 of the second parabola is such that both parabolas 
intersect at in the range < g, < 92- 



the ranges < q'l < qi and q2 < ^'2 < +00. Let us consider this probability in two steps. First, as for 
Ea. (|55p . we consider the initial conditions such that ipo{q) stays everywhere above a parabola Vxi,ci 
but goes below Vxi.ci+dci somewhere in the range < q'l < Qi- Integrating over the height ci this will 
take care of the first constraint < q'l < qi for the Lagrangian coordinate associated with xi. Second, 
we must only count among those initial conditions the ones that also satisfy q'2 > 92- We split them 
into two contributions as follows. Let us note q^ the unique abscissa where the two contact parabolas 
'Pxi,ci and Vx2,c2 intersect. From Eq.(fT2| it is given by 

X1+X2 C2-C1 
q* ^ — ^ 1. (81) 

2 X2 — Xi 

Then, we note p'' the first contribution, associated with initial conditions such that > (72 (which 
implies > 9* > 92)- Clearly, this actually corresponds to curves "ipoiq) that at some point go below 
the parabola Vx2,c, where c* is such that — q2 (i.e. the second parabola intersects 'Pa;i,ci at 92). We 
note the second contribution, associated with qi < q^ < q2 (since afterwards we shall consider the 
probability density Pa;i,x2('Zi; 92 — I2) we only need to include the cases with 9* > qi). We show in 
Figs, m and [5] the geometrical interpretation of these two contributions and p"^ . 

We describe in appendix [Cl the computation of the two-point distribution Pxi,x2{(li^ I2) from these 
two contributions and p^ . As for the one-point distribution computed in section 14.11 we first 
express the kernels K in terms of the Brownian propagator G obtained in section [3] and we use various 
properties of the Airy functions, described in Appendices 1X1 and [Bl to simplify the integrals. We finally 
obtain for the sum of both contributions the probability density 



Xl,X2\ 



+^°° d£ldS2 

(27ri) 



(Qi,Q2)=/ r:irrie(^i-i)'3i+(^^-i)'3^ij(si,2Xi)e-(v^~i)2X2i^ (82) 



Comparing with the one-point probability density (|73)) . we find that for < Qi < Q2 the two-point 
probability density factorizes as 

Xi<X2, 0<Qi<Q2: Px,^X2iQi,Q2) ^ PxAQi}Px2AQ2i), (83) 
where we introduced 

^e(-i)Qe-(yi-i)2X^ (84) 
-ioc 27ri 

Therefore, the conditional probability P{X2, Q2\Xi, Qi) obeys the property 

Xi<X2, 0<Qi<Q2: P{X2,Q2\X^,Qi)^^^^^^^^^§^^^Px2AQ2i), (85) 

PXi Wl) 
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that is, it only depends on the relative distances X21 and Q21, and no longer on Qi and Xi, over the 
range Qi > 0. Thus, the system is statistically homogeneous with respect to Lagrangian and velocity 
increments as long as we remain on either side of the origin Q = 0. Indeed, by symmetry through 
reflection about the origin we also have 

Xi<X2, Qi<Q2<0- Px,,x,{Qi^Q2) = P-x,,-xA-Q2,-Qi) 

= P-xA~Q2)Px,AQ2l) - Px,iQ2)Px,AQ2i). (86) 

Then, in the limits Xi —>■ +00, or X2 — > —00, where the weight of configurations such that Qi 
and Q2 have different signs vanishes, we recover the invariance through translation of the probability 
distributions of relative displacements and velocity increments. Of course, this is related to the fact 
that the initial conditions have homogeneous velocity increments, see Eq.®. However, the invariance 
through translation is only recovered in the exact nonlinear velocity distribution if we go infinitely 
far from the origin Q = 0. As we get closer to the origin it is broken by the increasing weight of 
configurations, such that Qi and Q2 are on different sides of the origin, which do not satisfy the 
factorizations ([55]) or ([55]) . Indeed, note that Eg. ([5^ shows that the factorization ([55]) cannot be 
extended to Qi < 0, as for Qi < and Q2 < we must reach the other regime (|86l) that cannot hold 
simultaneously. 

Thus, at finite distance from the origin the invariance through translation is always partly broken 
for the distribution Pxi,x2{qit Q2) considered over the full range —00 < 91 < 92 < 00. Nevertheless, 
the invariance is exactly recovered over either the partial range < qi < 172 < 00, or — 00 < gi < 
52 < 0. This can be understood as follows, focussing on the case gi > with again xi < X2- The 
probability density (qi) counts the configurations ipo{<l) that are tangent at qi with the highest first 
contact parabola 'Pxi,^, from the geometric construction described below (|12p . Then, the conditional 
probability density p{x2i,q2i\xi,qi) only counts among those the configurations that are also tangent 
at q2 with the highest first contact parabola Vx2,c2- Independently of the behavior of the curve ipo on 
either side of gi, the first-contact height parameter C2 must be smaller than the value c* such that 
'Px2,c, runs through the point {91, '0o(9i)}- Then, for any C2 < c*, we clearly have 'Px2.c2{q) < 'Pxi,ci{q) 
for all q < qi (using xi < X2), whence 'Px2,c2{l) < V'o(9) for all q < qi since we have already selected 
those configurations associated with Pxiiqi) that are above Vxi.a (and make contact at qi). (In other 
words, the additional requirement q{x2) = 92 does not bring any additional constraint on V'o('z) over 
q < qi.) Therefore, we are only sensitive to the behavior of ipo to the right of qi. For the Brownian 
initial conditions (O, the latter is fully determined by {?/'o(qi), i'o('?i)} and the white noise ^(g) at 
q > qi (which is statistically homogeneous). Next, p(x2i,q2i\xi,qi) does not depend on ipo{qi) since 
a vertical translation of the curves ipo and Vxi.ci is fully absorbed by the same vertical translation 
of the parabola 'Pxaxa: without affecting spatial coordinates q and x. On the other hand, through 
Galilean invariance the relative displacements of the particles only depend on their relative velocities, 
hence p{x2i, q2i\xi, qi) only depends on the relative velocity field VQ{q) — wo('Zi) over q > qi. For the 
Brownian initial conditions (O, with homogeneous velocity increments, the statistical properties of 
this relative velocity field wo(<z) — wo(9i) do not depend on vo(qi), but only on the distance q — qi, 
see ([9|). Therefore, the distributions p(x2i, 92112^1, gi), of the Lagrangian position increment (721, and 
p{x21tV2i\xi,vi), of the velocity increment V21, only depend on 0:21, as in ((85|) and in ([89|) below. 

In agreement with ([86)) . we can check that this argument fails for qi < 0. Indeed, again we are 
only sensitive to the behavior of iPo{q) to the right of gi, but this range now includes the special point 
q = with the constraints V'o(O) = and uo(0) = that prevent us from absorbing ipoili) and vo{qi). 
For instance, we now have the new constraint that the first-contact parabola Vx2.c2 cannot go upward 
of the point {O,-0o(O) — 0} (which was irrelevant in the previous case qi > 0, since we already had 
T^X2,C2 < 'Pxi,ci over q < qi and 7^a;i,ci(0) < by construction, being everywhere below ipo). 

The property that the increments of the inverse Lagrangian map, q{x2) — q^xi), are independent 
and homogeneous, as in Eq. ([85|l . and the probability distribution ([54|) . were already obtained by [13] 
for intrinsic statistical solutions, and by Q through probabilistic tools for a: > in the case of one-sided 
Brownian initial conditions (i.e. vo{q) = for q < 0). The latter work involves a similar reasoning to the 
one described above, using the property that the distribution of a Markov process after last passage at 
a given point does not depend on its previous path, but this mathematical proof uses the convex hull 
of the Lagrangian potential rather than the parabolas construction used here. For one-sided Brownian 
initial conditions, it is clear that if we consider Eulerian locations at a: > 0, the particles can only come 
from the right side q > so that we recover the configuration analyzed above for particles that are all 
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located on the same side of the origin. The agreement with the results of provides a nice check of 
our calculations. The probabilistic proof is remarkably concise, as it first shows that the increments 921 
are independent and homogeneous and next derives their distribution. However, the analysis method 
presented in the present work has the advantage of a large range of applicability. Thus, it allowed us to 
obtain the one-point distributions in closed form in section [4] and it could also be applied to different- 
time statistics, where the parabolas would have different curvatures. Another application of the method 
described in this paper is presented in [i^ , where we study ballistic aggregation for one-sided Brownian 
initial velocity. 

The previous discussion can be extended to n— point distributions, which thus factorize as 

Xi<..< x„, < Qi < .. < g„ : Pxu..,xAQu Qn) = PxAQi)Px,-xAQ2 - Qi) 

xPx.-xAQs - Q2)...Px„-X„_i(Q„ - Qn~l). (87) 



We obtain a similar identity for Qi < .. < Qn < by reflection through the origin, as for Eq. (|86|) . 
This also extends to the general case where the Lagrangian coordinates are located on both sides of 
the origin as 

X'^<..<X[<Xi<..< X„, Q'^<..<Q'i<0<Qi<..< Qn : 

m n 
1=2 j=2 

where we defined relative distances such as Xj^j^i — Xj — Xj^i. However, it appears that the probabil- 
ity distribution Px[,XiiQii Qi)i with Qi < < Qi, does not greatly simplify and is given by intricate 
multiple integrals. Therefore, we shall not consider it further in this article. Note that for practical 
purposes one is mostly interested in the behavior far from the origin, where the invariance through 
translations is fully restored. 

In terms of velocities, using the relation ([50)1 for the dimensionless velocities Vi, we obtain from the 
previous results the factorization 

x:^<..<x[<Xi<..<Xn, v;>x[, K,_i>x:,_i, Vi<X,, < : 

m n 

Px-x, = Px'^^x, {vi, ^1) n^^:-,. ^Ui) n ^^...-1 (^^-.^-i)' (89) 

i=2 j=2 

where the various factors are the velocity probabilities that may be obtained from the Lagrangian 
Q— probability densities through ((50|) . As noticed above, the factorizations (j88|l - (|89|) also follow from 
the analysis of However, although this provides the conditional distribution Px2-Xi{Q2 ^ Qi) 
of the Lagrangian increment it does not give the distributions Px'.:,Xj{Q'i',Qj) or -Px;,Xi (V^', Vi) that 
appear in these n-point distributions. Nevertheless, in the limit where we are far from the origin, we 
only need the one-point distribution Pxi{Qi), which goes to the Gaussian ([79|) . as would also be the 
case for one-sided initial conditions, besides in that limit we are mostly interested in the distributions 
of relative increments. 



We can note that the the Burgers ec 
paper was also used in a recent article 



nation with Brownian initial velocity which we study in this 
20l | to discuss the concept of local homogeneity that is used 
in turbulence studies. Indeed, for systems which are not strictly homogeneous (the energy shows an 
infrared divergence) it is customary to assume incremental homogeneity so that the physical quantities 
of interest (e.g. velocity increments) remain homogeneous. However, as noticed in 20] this is not fully 
consistent because initial incremental homogeneity is destroyed at later times by the nonlinearity of the 
equations of hydrodynamics (the quadratic advective term). Then, they used numerical simulations 
and perturbative analysis of the 1-D Burgers dynamics with two-sided Brownian initial velocity to 
illustrate this point and to note that local homogeneity is only asymptotically recovered far from the 
reference point. The results ([55)1 and ([55]) above explicitly show how the incremental homogeneity is 
indeed destroyed at finite distance from the origin but asymptotically recovered at large distances. A 
peculiarity of this system is that at finite distance it is already exactly recovered over a partial range 
of velocities. In fact, for the case of one-sided initial conditions (^0(9) = for g < 0) the system is 
exactly homogeneous over a; > 0, as shown by the previous discussion and 
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Fig. 6 (color online) Left panel: The probability density p^{q) that two particles, separated by the distance 
a; > at time t, were initially separated by a distance q (in the limit where the particles are far from the origin, 
or anywhere on the right side for one-sided initial conditions). This is the distribution of the increments of 
the inverse Lagrangian map, a; i— > q. We show the reduced probabilities, Px{Q), in terms of the dimensionless 
variables X — x/{2Dt^) and Q = q/{2Dt^), for three values of X, from Ea. (|90)) . The probability is zero for 
Q < 0. For large relative distance X we recover a Gaussian of center X and variance {{Q — X)^) — X/2. Right 
panel: The probability density Px{Q) on a logarithmic scale, for three values of X. 



The factorizations (|88p and ([55]) also show that small scales are largely decoupled from long- 
wavelength modes. Note that this key property is usually assumed in hydrodynamical systems (so that 
one can ignore the details of the large-scale boundary conditions) but is often difficult to prove in a 
precise manner. 



5.2 Distribution of Lagrangian increments (i.e. of relative initial Lagrangian distance) 

We now study in more details the probability distribution, Px{Q), of the relative Lagrangian positions 
(i.e. relative initial distance q between particles that are separated by distance x at time t), that is, 
of the increments of the inverse Lagrangian map (here we omit the subscripts "21" to simplify the 
notations) . The following results apply far from the origin, or at any location on the right side of the 
origin if we have one-sided initial conditions {vo{q) = at g < 0). 

We can check from the integral representation ((84|) that for X ^ we obtain as expected the Dirac 
distribution Px{Q) 5{Q) (whence Q2 Qi for X2 Xi). In fact, Eq. ([M)l is a well-known inverse 
Laplace transform [l| which gives the explicit expression 

X>0, Q>0: Px(Q) = AQ-3/2g2X-Q-XVQ = ^Q-3/2p-(VQ-X/VQ)^ (QQ) 

Therefore, we obtain an exponential tail at large Q, as ~ e~'^ , and a strong falloff at small Q, as 
~ e"'^ Z*^. For large relative distance X this gives the Gaussian 

\Q-X\^X: PxiQ)^ ^=e'^Q'''')"/''. (91) 

V ttX 

This agrees with the expectation that over large distances particles are still governed by the initial 
velocity field, as discussed in section [4?3l for Eg. ([79]) . This is again an illustration of the "principle of 
permanence of large eddies" ^3] : see the discussion below Eq. ([75)1 . 

We show the probability density PxiQ) obtained for three relative distances X in Fig. [51 We clearly 
see that for large X, which corresponds to large scales or small times, we recover a Gaussian centered 
on X, whereas for small X we obtain a skewed distribution with an intermediate power-law regime 
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From Eq. (190]) we obtain the moments of the Lagrangian increments Q as 22] 

^ _^ ^„«„ ^ g + ,92, 

where the last equality only holds for n > 1, and Ki, is the modified Bessel function of the second kind. 
This gives for the first few moments 

(Q>=X, (Q2)=X2 + |, (g3)=x^ + ^ + ^. (93) 

We can note that the mean of the relative displacement, x = X — Q, is zero: the mean distance 
between particles does not change (far from the origin). On the other hand, if we define the usual 
moment-generating function ^{y) by 



'^(y) = EM^(Q")= / dQe-^QPx(O), Px[Q)= ^cQy^iy) 
^0 Jo J-^oo 27ri 



(94) 



we obtain from Eq. (|84p . making the change of variable s = 1 + y, 

.Z/(y)=e-(^"i)2^. (95) 
Therefore, the cumulant-generating function ^(y), which satisfies the standard relation 

^(y) = E^(Q")c-M«'(2/)], (96) 

is given by 

<!>{y) = -{^TT~y l)2X ^ -Xy + 2X £ (-y)". (97) 



This yields the simple results 



(Q), = X, and for n > 2 : (Q"), = ^ X. (98) 

We can note that the first equality in ([5^ also holds for non- integer n, and we obtain for small 
Eulerian distance, {x2 — xi) — + 0+, 

u>\: {[q2-q^Y)^[2De)i-^)Et—^{x2-x^), (99) 
2 V"" 

^^<i: ((92 - qiT) - (2iJt^)-- ^^'"^ {x^ - Xi)^-. (100) 

Note that the second scaling also holds for any negative v. Indeed, the strong cutoff, e~'^ ^ of the 
probability distribution ([50]) . ensures that all negative moments are finite. Equations ([M| - (|100p show 
that we recover the bifractality of the inverse Lagrangian map, that was already derived in Q for v >Q. 
As is well-known [To^, the scaling ([M)) is universal as it is due to shocks. Indeed, if we have a shock of 
finite Lagrangian increment 5q at position a;, it gives a contribution [q{x -f t/2) — q{x — £/2)]^ ~ (^9)" 
which remains of order unity for £ — > 0+ for any n. Next, the probability to have a shock of a given 
finite strength 6q in a small Eulerian interval £ scales as £ at small distances, which gives rise to the 
factor (x2 ~ ^1 ) in Eq. ((99| . Note that in our case, the total number of shocks per unit length is actually 
infinite [4l|, |43|, see sect. IS.ll below. as the shock mass function (|154p leads to a divergence at low mass, 
but the number of shocks above a finite mass threshold is finite and this is sufficient to make the scaling 
([M)) universal. However, the behavior observed at < 1 (the critical value i^c — 1/2 and the exponent 
21/ observed below Vc in Eq. (jlOOp ) depends on the initial energy spectrum, through the low-mass tail 
of the shock mass function, see also ^] for more detailed discussions. 
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Fig. 7 (color online) Left panel: The probability density p^{v) of the velocity increment v — v{x2) — v{xi) for 
two positions separated by the distance x = X2 — Xi (in the limit where we are far from the origin, or anywhere 
on the right side for one-sided initial conditions). We show the reduced probabilities, Px{V) in terms of the 
dimensionless variables X = x/{2Dt^) and V — v/{2Dt), for three values of X, from Ea. (|101|l . The probability 
is zero for V > X. At large relative distance X we recover a symmetric Gaussian of variance 
Right panel: The probability density Px{V) on a semi-logarithmic scale, for three values of X. 



5.3 Distribution of Eulerian velocity increments 

We now consider the probability distribution, Px{V), of the relative Eulerian velocities, that is of the 
velocity increments V{X2) — V{Xi). From Eq. ([90l) we obtain 

X>0, V<X: PxiV) = ^iX-V)-^/^e-^^^-''/^^^\ (101) 

In the limit of large relative Eulerian distance X — > oo, at fixed V, the distribution (|10ip can be 
expanded around the maximum of the exponent at = (corresponding to Q = X) and we again 
recover the initial Gaussian 

\V\^X: Px(l/) ^ _^e-^'/^, (102) 
V ttX 

in agreement with the fact that over large distances particles are still governed by the initial velocity 
field (see also section We show in Fig.[7]the velocity distribution PxiV) for three values of X. We 
can again check that we recover a Gaussian for large X (i.e. large scales or small times), whereas for 
smaller X the upper bound V < X is increasingly apparent while a power law develops at intermediate 
negative velocities. 

From Eg. ([55)1 the velocity cumulants are given by 

n>2: {V")c = (-l)" ^^^~f' X, whence {V) = 0, (V) ^ |, {V') = (103) 

We can note that the first moment exactly vanishes whereas the variance {V'^) remains equal to that 
of the initial Gaussian field, see ([HI), even though Px{V) is no longer Gaussian. Thus, in terms of the 
dimensional variables the velocity energy spectrum remains equal to the initial one, 

{[v{x2,t)-vixi,t)]^) ^D\x2--xi\, E{k,t) ^ Eo{k) ^ ^k-\ (104) 

Finally, in the limit of small separations we obtain from Eq. (|103[) 

n>2, (x2-x^)^0+: {{v^ ~ ViT) - ^ J (-1)" ^ (^2 - xi). (105) 

Therefore, we recover the universal scaling at small distances of the structure functions jT^, {[v{x + 
€} — u(x)]") cx £, that was also observed in the numerical simulations of |41|]. This is due to the 
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Fig. 8 (color online) Left panel: The probability distribution Pxirj) of the overdensity, r] — m/{pox), over a 
region of length x. We show Px{ri) for three values of the reduced length X = x/{2Dt'^), from Eq. p07|) . Thus 
larger X corresponds to larger scale or smaller time. For large X we recover a Gaussian of mean 1 and variance 
((»? - 1)^) = 1/(2^)- For small X the distribution becomes skewed and an intermediate power-law region 
develops. Right panel: The probability density Pxiv) on a logarithmic scale, for three values of X. 



contribution from shocks, as discussed below Eqs. ([99|) - (|100|) . Thus, if we have a shock of finite velocity 
jump -6v = Sq/t at location x, then [v{x + i/2) - v{x - £/2)]" - {-Sv)"^ for £ 0+. Note that Sv 
is positive, since a shock is associated with particles from the left overtaking particles on the right, so 
that v{x~) > w(a;+), which agrees with the factor (—1)" in Eq. (|105p . Again, the factor {x2 — Xi) in 
Eg. ljlOSp comes from the probability to encounter a shock of strength larger than some finite threshold 
Sq in a small Eulerian interval [xi, X2]- 



6 Density field 

6.1 Overdensity within finite size domains 

We consider here the evolution of a density field p{x, t) that evolves through the usual continuity 
equation, 

whereas the velocity field v{x,t) evolves through the Burgers equation ([T]). The initial conditions for 
the velocity are set by Eqs.® as in previous sections, whereas the initial density is a constant p^. Thus, 
the mass m between particles qi and (72, with qi < q2, is m — pn{q2 ^ qi)- This quantity is conserved 
by the dynamics since particles do not cross each other (though it is ambiguous at shock locations, 
but the latter have zero measure in Eulerian space). Then, the overall overdensity, 77 = m/{pox), over 
the length x — X2 ^ xi, is rj = ((72 ~ 'li)/{x2 ~ x\) by conservation of matter, where qi is the initial 
Lagrangian position of the particle that is located at Xi at time t. In terms of dimensionless variables 
this reads as the ratio of relative distances 77 = Q/X. Therefore, far from the origin (|a;i| 00), or 
on the right side of the origin for one-sided initial conditions, we obtain from Eq. (|90p the probability 
distribution of the overdensity at scale X as 

^=^, ^>0: P^(^)^./Z^-3/2e-x(v^-i/v^)^^ /^g2X^-3/2g-x(,+iA,)^ (^07) 
Pqx V tt V tt 

Over large scales we recover a Gaussian distribution, as for the variables Q and V , while on small scales, 
X — > 0, we obtain the power law r\~'^l'^ between the low and high density cutoffs at r\_ ^ xj (2Dt^^ and 
77+ ~ {2Dt^)/x, as we can see in Fig. |8] where we show the overdensity distribution Pxiv) for three 
values of X. We can note that this is very similar to the behavior that is observed in cosmological 
numerical simulations for gravitational clustering 0, Hi 113 ■ 
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From Eqs. ([92|) and ([98|) we obtain for the moments and the cumulants of the overdensity at scale 

X: 

whence for the lowest orders 

(^> = 1, = (^'>c = ^ = 3(^2)^ (109) 

We can note that the second result (|108p gives the cumulant hierarchy 

Sn = fy^^{2n-i)\\ and ^(y) ^ ^„ ^ = yrT2^ - 1, (110) 

which shows that the ratios Sn are constants that do not depend on time nor scale. We can note that 
in the cosmological context, associated with a gravitational dynamics in a 3-dimensional expanding 
Universe, for the case of an initial power-law density power spectrum, the coefficients 5*^, still defined 
as in (jllO[) . only show a weak dependence on scale in the highly nonlinear regime, and they also 
asymptotically reach (different) finite values at large scales in the quasi-linear regime [11,11. Thenjt 
has been proposed to use the approximation of constant Sn to describe the highly nonlinear regime [37| . 
Moreover, the form ()110|) of the reduced cumulant generating function (/7(y) is one of the possibilities 
that have been studied in this context This phenomenological ansatz is known as the "stable 
clustering model" as it was derived by assuming that on small physical scales, after nonlinear collapse 
and gravitational relaxation, overdensities decouple from the Hubble expansion and keep a constant 
physical size [l5|. In the present case, collapsed objects are actually Dirac peaks (shocks) of vanishing 
size. Then, it is easy to see from a multifractal analysis that shocks lead to finite ratios Sn in the small- 
scale limit llll, so that the hierarchy (|110[) is universal, in the sense that the generating function 
^p{y) has a finite limit at small scale, x Q. However, this non-trivial limit depends on the initial 
energy spectrum. A specific property of the Brownian initial conditions studied in this paper is that 
the ratios Sn are actually constant over all scales, from the linear to highly nonlinear scales. Thus, it 
is interesting to note that the 1-D Burgers equation with Brownian initial velocity provides an exact 
physical realization of this ansatz. 



6.2 Density correlations 

We now consider the unsmoothed density field p{x) itself (again in the limit where we are far from 
the origin so that the previous results apply). It is related to the smoothed overdensity rj over scale x 
introduced above through 

X pa 

Then, introducing the density power spectrum, V{k), by going to Fourier space as 

p{x)^Po= f dfce"=^p(fc), {p{ki)p{k2)) ^6{k,+k2)plr{k,), (112) 

J —oc 

we obtain using the second Eq. (|109p and Eas. (|llip - pi2p 

Df^ r°° kr Bf'^ 

— ^{v')c^ dksmc'{—)V{k), and V{k,t) ^ —, (113) 
X 2 27r 

where sinc(a;) — sin(a;)/x is the cardinal sine. Thus, we obtain a white-noise density power spectrum, 
with an amplitude that grows as t'^. This yields the connected density two-point correlation 

{p{xi,t)p{x2,t))c ^ plC2ixi,X2), with C2{xi,X2) = Dt^ d{x2 - Xi), (114) 

which remains a Dirac function at all times. 
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+ 5 terms + (2 J (4') + 2 terms + (4j (ij + 2 terms + ^2) + 1 term + 

3 ) 13, 



Fig. 9 The 15 heap ordered trees that can be associated with the 4— point correlation C4,{xi, .., X4). The labels 
refer to the positions Xi, which are ordered as xi < X2 < X3 < Xi. We only show the 5 tree structures, as the 
additional terms can be obtained from the previous diagram by permutations over the labels that satisfy the 
ordering constraint that each path from the root has increasing labels as we proceed down to the leaves. Each 
link between nodes i and j yields a contribution C2{xi,Xj) and the contribution of a tree is the product of the 
3 factors C2 associated with its 3 links. This gives C4 as the sum over all these tree contributions. 



In fact, the factorization (|87p implies a similar factorization for the multivariate distributions of 
the smoothed density field, far from the origin (Xi +00), 

Xi < .. < Xn : Px2,i;..;X„,„„i(?72,i; •■;'7n,»i-l) = -Pxa.i (^2,l)--^'x„,„_i (f?n,«-l), (115) 

where 77i,i-i is the mean overdensity over the interval Thus, the densities within non- 

overlapping domains are completely independent random variables. This agrees with the Dirac obtained 
in Eq. pi4p for the connected density two-point correlation. Moreover, this can be extended to all higher 
orders. Indeed, let us consider the density n— point connected correlation, defined as 

{p{xi)..p{Xn))c = Po C'nixi, ..,X„). (116) 

If there exists a position Xi that is different from all other positions Xj, with j ^ i, then we can build 
a small region [xi — e,Xi + e] where the density is independent from the density at all other points 
Xj, using the property pi5|) and e — > 0"*". Therefore, by definition of connected correlations, C„ must 
vanish. Then, the point connected correlation can be written as the product of n — 1 Dirac factors 

Cn{xi, ..,Xn) - {2n - 3)!! (Dt^)^-^ Six2 - Xi)6{x3 - xi)...5(:e„ - xi) (117) 
= {2n ^ 3)l\C2{xi,X2)C2ixi,X3)...C2{xi,Xn), (118) 

where the amplitude is obtained from Eq. JlOSp . as well as Eq. piip which implies the relation (77") c = 
x~" dxi..da;„C„( Xi , . . , Xn } . 

We can note that (2n— 3)!! also counts the number of heap ordered trees with n nodes, that is, rooted 
trees where the n nodes are labelled as {1,2, ..,rt} and each path from the root has increasing labels 
[39| . In our case, we can therefore construct the following combinatorial interpretation of Eq. (|117p . 
First, the points a;i,..,x„, are ordered such that < X2 < .. < Xn (we choose one among several 
possibilities if several positions are equal). Then, the n— point connected correlation ()117|) is obtained 
as the sum of the contributions of all heap ordered trees, where the contribution of each tree is simply 
the product of the n — 1 factors C2{xi,Xj) associated with the n — I links between nodes Xi and Xj. 

Of course, we may also write (I117P as the sum over the products of C2{xi,Xj) associated with 
any other class of N trees, multiplied by a weight {2n — 3)!!/A^. However, this no longer recovers the 
amplitude (2n — 3)!! of Eq. (|117p in a natural manner. We note that in the cosmological context, within 
the stable-clustering ansatz discussed above, it has been proposed to use as a phenomenological model 
a diagrammatic description such as Fig. [21 where the connected n— point density correlation is written 
as the sum over trees of each product of {n — 1) factors C2{xi,Xj) associated with the internal links 
[UliOl. However, the tree diagrams used in this context are usually not ordered, and each topology may 
have an additional multiplicative weight. In the present case of the 1-dimensional Burgers dynamics, 
we can note that the concept of ordering naturally arises since particles do not cross each other and 
one can order both Lagrangian and Eulerian positions on the line (this would no longer be the case 
for higher dimensions). 

Thus, the 1-D Burgers dynamics with Brownian initial velocity provides a physical realization of 
the hierarchical structure such as (|118p for the many-body correlation functions. It is an interesting 
question to ask whether other real dynamical systems can be built that display the same factorization 
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property (possibly over other classes of trees) with other two-point correlations C2EI On the other 
hand, one may wonder whether a factorization such as (jllSp . and a diagrammatic construction such as 
Fig. [9l could be generalized, as an exact asymptotic solution or as a useful phenomenological model, 
to the 1-dimensional Burgers dynamics with other initial conditions, where C2 would no longer be the 
simple Dirac function (jll4p FI 



6.3 Comparison with a perturbative approach 

We can note that the exact (far from the origin) nonlinear results pi3p - ()114p are identical to the 
perturbative predictions that would be obtained at linear order from Eq. (jl06p . Indeed, if we linearize 
the continuity and inviscid Burgers equations, we obtain at lowest order for the density field dpL/dt = 
—podvL/dx — —podvo/dx — — poCi where ^ is the initial white-noise of Eq.([5]). This gives pL{x,t) — 
po{l — t^{x)), which leads in turn to Eqs. (|113p - (|114p . The fact that for these Brownian initial conditions 
the nonlinear Burgers dynamics preserves the linear density power spectrum is reminiscent of the 
invariance of the energy velocity spectrum (|104p . In both cases, one needs to consider higher-order 
correlations (or the full distribution) to measure the effects of the nonlinearities. 

In fact, the agreement of the exact density two-point function with the linear theory actually 
extends to all order cumulants {rf')c, computed at leading order from quasi-linear theory. Indeed, at 
tree-order in perturbation theory, in the inviscid limit, it can be shown that the cumulant-generating 
function <p(j/), defined as in Eq. (jllOp . is given by the implicit system 



T -vQ'{t) 

with Q{t) = T 



c{G x) 

-T- 



a{x) 



^1/2 



(119) 



where (t(x)^ = {5\) = Dt^ /x is the variance of the linear density contrast Sl at scale x, and the function 
JF((5l) describes the evolution of spherical (here symmetric) density fluctuations (see 0,11,111] for the 
similar case of the cosmological gravitational dynamics). The system {r, Q} <-> {y, ip} in (|119p is actually 
a Legendre transform and it arises from a saddle-point approximation. Indeed, in the quasi-linear limit 
(i.e. tr — *■ 0, which also corresponds to t — > or x ^ 00) the cumulant ratios Sn are governed by the tails 
of the density distribution and the generating function ip{y) can be obtained from a steepest-descent 
methocfl [45!|. (In a somewhat similar fashion, the minimization problem ([3]), that also arises from a 
saddle-point method, can be written in terms of a Legendre transform of the Lagrangian potential, 
see p.) As compared with Eq.(69) of Q we made the change a[{l + G^^^x] a{Q x) by taking 
1 + Q ^ Q and the exponent 1/3 is changed to unity as we go from 3-D to 1-D. For the present 1-D 
Burgers dynamics, we have: 



■!]= — = — - — , whence at linear order rj^ = 1 and Sl = —t—. (120) 

X q + tv q q 



This yields 



^(^l)--^, whence g(r) ^ ^ ^ + + 4)' ^{y) ^ ^lT2^ - 1. (121) 

1 — ol 4 

^ In fact, as for the weaker property of constant ratios Sn, the author is not aware of other dynamical systems 
that exactly obey such a factorization property. In view of the many phenomenological works that have used 
such a diagrammatic construction for many-body correlations, it is satisfying to find that it is at least obeyed 
by one truly dynamical system, even though the expression in terms of Dirac factors and the lack of large 
distance correlation make this a very simple and specific case. 

The results of 9] show that the Brownian case can be generalized to Levy processes with no positive jumps, 
where the increments of the inverse Lagrangian map again remain homogeneous and independent at all times. 
However, this does not provide another hierarchy for the many-body correlations as they remain of the Dirac 
type. 

^ In fact, the steepest-descent method described in 0^ is a non-perturbative approach, which can als o be 
applied to the other limit of rare events at finite a, where it allows to go beyond perturbative methods [4g |. 
However, in the quasi-linear limit ct ^ it gives the same results for 'fi{y) as the usual perturbative expansion 
over powers of the linear growing mode of the density field (provided the latter gives finite results in this limit). 



24 



Thus, we recover at tree-order the exact resuh ()110p . Therefore, for Brownian initial velocity the 
Burgers dynamics happens to preserve the density cumulant-generating function (p{y) that would be 
obtained at leading order (tree-order) from a perturbative approach, which does not take into account 
collisions between particles (shell-crossings). This is also why the ratios Sn are constants that apply 
to all scales, from the quasi-linear to the highly nonlinear scales. Note that the perturbative approach 
breaks down beyond leading order as next-to-leading corrections actually involve divergent integrals 
(which means that one can no longer discard shocks, which requires non-perturbative methods). For 
other initial conditions the coefficients Sn would no longer remain equal to their tree-order values. 
However, they still asymptote to finite values in the highly nonlinear regime, because of the contribution 
from shocks, just as the Lagrangian and velocity increments scale linearly with £ for small distances 
^ ^ 0, as discussed below Eas. ((M)) -([TOO | and Eg. pOSl) . 



7 Lagrangian displacement field 

7.1 One-point distributions 



We now consider the dynamics associated with the Burgers equation H]) from a Lagrangian point of 
view. That is, labelling particles by their initial position q at time t = 0, we follow their trajectory 
x(q,t) and we note xilit) — ^{lit) ~ Q their displacement with respect to their initial location. Note 
that for regular points, which have kept their initial velocity, we have x = see Eq.Q. Since particles 
do not cross each other they remain well-ordered. Then, it is clear that the probability, Pq{x' > x), 
for the particle q to be to the right of the Eulerian position x, at time t, is equal to the probability, 
Px{q' < q), for the Eulerian location x to be "occupied" by particles that were initially to the left of 
particle q. (Since shocks have zero measure in Eulerian space there are no ambiguities.) Therefore, we 
obtain in terms of dimcnsionless variables, for the case q > 0, 



g > : PQiX'>X) = PxiQ' <Q) = PxiQ' < 0) + Px(0 < Q' < Q) 



/+ioo 
--ioo 



ds_ 
2ni 



J{s,2X) 
s-1 



(122) 



where we used the results of section UT] and the integration contour runs to the right of the pole s = 1. 
Therefore, the probability density of the Eulerian position X of particle Q reads as 



Q>0: Pq{X) 



d 



dX 
d_ 

' ^dX 



Pq{X' > X) 



J{l,-2X)+ I f:e(-i 
toe 27ri 



)Q J{s,2X) 
s- 1 



(123) 



The case Q < can be obtained from Ea (|123p through a reflection about the origin. 

Let us consider more precisely the case of the particle that was initially at rest at the origin, 
q = 0. Following the previous discussion, we have pq{x' > x) — Px{q' < 0), hence we can directly use 
Eqs.(I7Dl)-(172]) which give 



X>0: Po{X) 



d 



dv Ai 



dX Jo V 

At small X, we obtain from Eq. (|124p the asymptotic 



v2X 



Po{-X) - Po{X). 



X^0+: Po{X) ^ -—\nX, 



(124) 



(125) 



whereas the behavior for large displacement X is set by the asymptotic ([50)) . which yields 



X 



+00 : PoiX) - (7rA:/(2%/3)) e 



(126) 



Thus, the central particle q ~ 0, that was initially at rest {vq{Q) — 0), has moved by time t by a distance 
X = X whose distribution shows an exponential tail at large \X\ and a logarithmic peak at low \X\. 
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X = x/(2Dt) log[X]=log[x/(2Dt)] 



Fig. 10 (color online) Left panel: The probability distribution, Pq(X ), of the reduced position X = x/{2Dt) 
of the par ticle that was i nitially at the origin, q = 0, from Ea. (|124|l . The dashed lines show the asymptotic 
behaviors (|125p and (112611 . Right panel: Same as left panel but on a logarithmic scale. 



We can note that both the large-X tail and the low-X divergence are different from the asymptotics of 
the distribution of the Lagrangian coordinate, Q — — V , oi the particle located at the Eulerian location 
X = at the same time, see Eqs. (l75)) - ([77)) . We show in Fig. [TU]the distribution Po{X) as well as the 
asymptotic behaviors (|125p and (|126p . It appears that the logarithmic asymptote is only reached at 
very low X. 

Finally, far away from the origin, in the limit Q — > oo at fixed x = ^ ~ we obtain from Ea. (|123p 
the asymptotic behavior (making the change of variable s = 1 + ifc as for Eq. ([79| ) 

Q^+oo, |x|«Q: Pq{x)-^^. (127) 

Therefore, we recover the property that far from the origin the displacement of the particles is governed 
at leading order by the Gaussian distribution of the initial velocity field vo{q)- This agrees with the 
discussion presented below Eg. 179]) in section l473l Again, this expresses the fact that at very large 
distances, where the initial velocity becomes increasingly large as y/\q\, the motion with respect to the 
origin is dominated by the "large-scale flow" and the local fluctuations of the initial velocity field have 
only produced local subdominant shifts, see also [24>. ,3.] . 



7.2 Two-point distributions 



We now investigate the two-point probability distribution of the Lagrangian displacement field. As for 
the one-point distribution (|122p . it is related to its Eulerian counterpart through 

PQr,QAX[ >X,,X'2> X2) - Px^,xM <Ql,Q2< Q2) (128) 

Using Eas. ([75)) . ([M]) . we obtain far away from the origin, in the limit Qi +00 at fixed Q21 — 
Q2 - Qi > 0, 



?2i>0: Pq,^qAX[>X,,X'2>X2) 
dk 



Qi 



Ak{Xi-Q[)-Xik^/4 



2tt 



+ZOO 



Q2 

dQ[ I dQ'2 
JQ'i 

ds 
27ri 



g(s-l)(Qi-Q'i) g-(V7-l)2X2i_ (^129) 



Integrating over Q'l and Q2 gives the cumulative distribution 

r Ah r+ioc 1 

Qi ^ Q21 > : Pq„qM'i>Xi,X'2>X2)^ — — 

J 27r J_i^ ln\ 



-ifcXi-Xifc^/4-(ys-l)2X2i r 



(s-l)(i/c 



1) 



„ifcQl + (s-l)Q21 _ „(s-l)Q2 



(130) 
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Then, taking the derivatives with respect to Xi and X2, and integrating over k, gives in the hmit 
Qi — !■ +00, at finite xi = ^1 ^ Qii the probabihty density 

VT^Qi J-joo 27ri (Vs + 1)^ 

Thus, the comparison with Ea. p27p shows that we obtain as expected a factorization of the form 

gi -> +00, Q21 > : Pq,,q,{Xi,X2) ^ PqAXi)Pq,AX2i), (132) 

where the distribution, Pq^-^{X2i), of the relative Eulerian distance X21 of the particles that were 
initially separated by the distance Q21 reads as 

.... S^"-""-'^-""'l7JTIF- ''''' 



Therefore, we recover a factorization of the form (|83p that was obtained in the Eulerian framework. 
However, it is no longer exact at a finite distance from the origin and only applies in the limit Qi — > 00 
(again, by symmetry we have a similar result for Q2 — > —00). 

Multiplying Ea. (|133p bye"'*'^ and taking the derivative with respect to X we obtain a standard 
inverse Laplace transform [ij 



1 <.+ioo 1 o -(v^+l)2X 



(s-i)Q 



= -8 



(134) 



where erfc(z) is the complementary error function. This can be integrated to give 



Pq{X) ^8{X + Q+^)e^'' erfc ( ^ + j - 8 y|e"<^^-^'^ (135) 



Using the asymptotic expansion of the complementary error function [ij we obtain for large Lagrangian 
separation, Q, and fixed relative displacement, % = X — Q, 

Pq{x)^ ^e'^"/Q. (136) 

As for Eqs. ((79)l . (|102p . (|127[) . we recover as expected the property that over large distances particles 
are still governed at leading order by the initial Gaussian velocity field. Next, equation (|135p yields for 
the asymptotic behavior at large X for finite Q, 

X^+00: PQ{X)^\[^-e'^V3-^^' =^f^erQx-'e^''~'''^Q, (137) 

\ TT X V TT 

whereas Pq{X) remains finite for X 0+. We show our results for three values of Q in Fig. [Til that 
clearly illustrate the evolution of Pq{X) with scale or time (smaller Q corresponds to smaller scale or 
larger time). As for the Eulerian probability distribution, the Gaussian tail (I137P can be understood 
from a simple scaling argument applied to the initial velocity field. Thus, the expansion of the initial 
Lagrangian interval g up to a very large size x at time t requires an initial velocity increment of order 
vo ~ x/t (since x ^ q) which gives rise to a probability of order e~^^/*^ ~ e~'^ Z*^, using Eq.®, 
which agrees with the large-X tail (|137p . 

From Eq. (|135p we also obtain for the asymptotic behaviors of Pq{0^) with respect to Q 

Q^O: Pq(0+)->4, +00 : Pq(0+) - ^g~3/2e-Q. (138) 
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Fig. 11 (color online) Left panel: The probability density Pq{x) that two particles, that were initially separated 
by a distance q, are separated by the distance a; > at time t (in the limit where the particles are far from the 
origin). We show the reduced probabilities, Pq{X), in terms of the dimensionless variables Q = q/{2Dt^) and 
X = x/{2Dt^), for three values of Q, from Ea. (|135|) . The probability is zero for X < 0. For large initial relative 
distance Q we recover a Gaussian of center Q and variance {{X — Q)^) = Q/2. Right panel: The probability 
density Pq(X) on a logarithmic scale, for three values of Q. 
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Fig. 12 (color online) Left panel: The probability, p"^"'^^, that two particles, that were initially separated by 
a distance q, have coalesced within a single shock by time t (in the limit where the particles are far from the 
origin). We show the reduced probability, Pq°'^'^, from Ea. (|144p . The dashed lines are the asymptotic behaviors 
(11421) and (|143|) . Right panel: The probability p*°'^'^ on a logarithmic scale. 



However, note that the limits X — > and Q — > do not commute. Indeed, it is clear from Eq. (|135p 
that for any X > we have Pq{X) for Q — > 0. As we shall see below, this is the signature of the 
contribution due to shocks. Thus, from Ea. (|133p we obtain the cumulative distribution as 



g > 0, X >0: Pq{X' > X) 



-\-ioo 



ds 
27ri 



(s-l)(Vi+l)' 



(s) > 1, 



where the integration path is located to the right of the pole at s = 1, so that Pq{X' > X) 
X —>■ +00. However, we note that it does not reach unity in the limit X 0^ , since we have 



(139) 
. for 



Q > : ^lim^ Pq{X' > X) = 1 



-\-ioo 



ds_ 
27ri 



, < 5R(s) < 1, 



(140) 



where the integration path crosses the real axis in the range < s < 1. This means that there is a 
non-zero contribution due to shocks, where particles that were initially located at different positions, 
qi ^ q2, have collided by time t and are now located at the same Eulerian position, xi = X2, in the 
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same massive shock. Therefore, to the contribution ()133|) we must add the contribution from shocks, 
that reads as 

g > : Pq°''^{X) = Pj""*" S{X - 0+), with the amphtude 

^e(-^)^ ^ 0<5R(.)<1, (141) 

-joo 27ri (l-s)(Vs + l) 

so that the full probability is normalized to unity, see Eq. (|140p . Thus, the amplitude p^°'^'* jg 
probability that two particles, that were initially separated by the (dimensionless) distance Q, are 
both located in the same shock at time t (in the limit where the particles are far from the origin, 
or anywhere on the right side for one-sided initial conditions). At large initial Lagrangian separation 
Q, we obtain from Eq. (jl4ip the exponential decay (that again can be understood from simple scaling 
arguments) 

Q^+oo: pj"^"^^ 1 Q-3/2g-Q^ (142) 
whereas for small initial distance Q we have 

pJ°^'^^l-4./:^. (143) 

V TT 

Therefore, in the limit Q — s- the probability that both particles are within the same shock reaches 
unity whereas the weight associated with the "regular" contribution (|135p vanishes (while its cutoff 
decreases as Q). This agrees with the well-known result that the set of regular Lagrangian points has 
a Hausdorff dimension equal to 1/2 [4l|, |4^, so that with probability 1 a random Lagrangian point 

q belongs to a shock at any given time t > 0. This clearly implies that Pq — > 1 for Q 0, as 

in Eq. (|143p . Moreover, the behavior 1 — p^^°'^^ Q^^^ also shows that the set of regular Lagrangian 
points has a box-counting dimension equal to 1/2, in agreement with these works. 

Taking the derivative of Eq. ()14ip yields again a standard inverse Laplace transform 1] that provides 

1 • r -prshock 

a convenient integral expression for Pq , 



—shock 



dO 



erfc(VQ) 




hence pT = 2 / dQ' -y= - erfc(\/Q^) . (144) 



We show in Fig. [T^the same-shock probability p^^""^^ a,s a function of Q, as well as the asymptotic 
behaviors (fT42| and pl3| . 



7.3 Higher-order distributions 

We can obtain the higher-order n— point distributions Pq-^ . . . ^g,^ (xi , . . , x„ ) by the same method which we 
applied in the previous section for the two-point distribution. Thus, as in Eq. (|128p . we can related the 
Lagrangian and Eulerian cumulative probabilities by 

PQ^,...QAX[ > > X„) = Px„...xAQ'l < Qu-,Q'n < Qn), (145) 

with Qi < Q2 < ■■ < Qn and Xi < X2 < .. < X„. Then, in the limit Qi +00, using again the 
factorization ([57)) and the expressions ([75)1 and ([M)l . we can integrate over Q[, Differentiating 
with respect to Xi, .., X„, gives the 71— point probability density (compare with Eq. ()13ip ) 



Q.^oo: Pg„..,g„(X„..,X„)-^=^y_ J^^^ ' 



on „-'(V^-l)2X2,i-..-(Vs^-l)2X„,„_i 

— —— , (146) 

(1 + ^/S^){^/S^ + y/s^)..{y^Sn-l + ^/S^){y/S^+ 1)' 
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with Xi^i-i — Xi — Qi^i-i ^ Qi — Qi-i, ^1 = Xi — Qi. Note that this n— point distribution does 

not factorize. 

From Eq. (|146p we can obtain the contributions associated with shocks in the same manner as in 

section [7^ For instance, far from the origin [Qi ^ 1), the probabihty density, -Pq2°%32,Q43 (^32)1 that 
each pair {Qi, (52}, and {Q3, Q4}, has coalesced within two shocks that are separated by a distance in 
the range [X^2,X^2 + dX32], reads as 

—shock /• + '°° dS2dS3dS4 (^2-l)Q21 + (^3-l)Q32 + (^4~-l)Q43 

4e-(v^-i)2^32 



(s2 - l)(ys2 + V^)('54 - 1)' 

where the integration contour is such that 3fJ(s2) < 1, 5R(s3) > 1 and '^{34) < 1. 



(147) 



7.4 Computing the density power spectrum from the Lagrangian statistics 

FinaUy, it is interesting to note that the statistics of the Lagrangian displacement field also allow us to 
compute the Eulerian density power spectrum and to recover the result (|113p . Indeed, as is well-known 
the conservation of matter implies that the density p{x) may be written as 

p{x) = Pa f ^) =po f dqdix-q-- x(g)), (148) 



where x{q) = ^il) — 9 is the Lagrangian displacement of particle q. Note that the last expression is 
still valid when there are shocks, as may be seen by computing the mass within some interval [a;i,a;2]. 
Going to Fourier space as in (|112p . we can write 

(p(fci)p(fc2)) ^pI [ ^§l^e-('=^^^+'=^^^)(e-'('=«^+'^^^^)). (149) 

Then, in the regime where the invariance through translations is recovered (i.e. far from the origin), 
making the changes of variables q2 = qi + q and X2 = Xi + X: we obtain 



f°° dq f°° 

/ — / d.Tp (a;) cos(fca;) (150) 
Ja Jo 

1' r- rdXPQiX)cos{^'kX)+^^ rdQ-s^o.. ^^^^^ 

Jo Jo+ Jo 



V{k) 



In the second line, written in terms of dimensionless variables, we separated the two contributions 
associated with the regular part p35p and with the singular part (|14ip of the distribution of the relative 
Eulerian distance X. Using Eq. (|133p . we can check that the first contribution actually vanishes (as 
expected since all the mass is enclosed within shocks, see Eq. (jl53p below and [4l|, [43!|), whereas the 
shock contribution is obviously independent of k and using Ea. (|14ip we recover the amplitude (|113p . 



8 Properties of shocks 

8.1 Shock mass function 

From the probability jf^°'^'^ that two particles, initially separated by a distance q, are located in the 
same shock at a time i > 0, we can now derive the mass function of shocks. First, we note that if a shock 
has Lagrangian end-points g_ and q+, its mass is simply m = po{q+ — q-), where po is the uniform 
initial density, as discussed in section [6l Then, within an interval of length Q in Lagrangian space, we 
now count the number, Quq (m)dm, of shock-intervals of length in the range [g, q + dg] , whence of mass 
in [m, m -|- dm], with m = poq. The limit Q — *■ 00 gives the probability density n{m) for a Lagrangian 
point to belong to a shock of mass m. Since on large scales particles are still governed by the initial 
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Gaussian velocity field and have only moved by a relative distance x ^ VQ, see for instance Eq. (|136p , 
the corresponding Eulerian relative distance, X, obeys X/Q ^ 1 for Q ^ oo. Therefore, n(m)dm is 
also the mean number of shocks, per unit Eulerian length, with a mass in the range [m, m + dm]. 

Let us now relate the mass function n{m) to the shock probability that we obtained in 

Eq. (|14H ). The latter is the probability that a Lagrangian interval, q = — qi, chosen at random, has 
coalesced by time t within a single shock. This clearly means that this shock has a length larger 
than q and that qi is located within a distance smaller than q^ — q from its left boundary. Therefore, 
in terms of dimensionless variables, we have the relation 



Po°'''= / dMN(M)(M~Q), with Af=-^. (152) 

PoT 



Here N{M) is the dimensionless mass function. From Ea. (|152p and Ea. (|143p . taking Q = 0, we obtain 
at once the normalization 

/•oo 

dMMN{M) = l, (153) 







which means that all the mass is included within shocks, at any time t > Q. This agrees with previous 
results discussed below Eq. p43|) . see [4l| and [43|. Differentiating twice Eq. p52|l with respect to Q, 
and using the first Eq. p44p . gives the simple expression 



N{M) - 



dg2 



whence N{M) ^ A-r^''^ e'^' . (154) 



Thus, we recover the low- mass power law M~^/^ that was already obtained in [4l|, [4^. At large 
masses we obtain the exponential falloff that was heuristically derived in [5l|, following the same 
scaling arguments as those described in previous sections for the tails of the Eulerian or Lagrangian 
distributions. The full mass function (|154p was also obtained in [§] for the one-sided Brownian initial 
velocity. Indeed, in that case the Lagrangian increments q2i have the same distribution for a; > 0, as 
discussed in section [5Tl which clearly leads to identical shock properties. 

It is interesting to compare the exact result (I154p with the Press-Schechter ansatz that is widely 
used in the cosmological context to count the number of collapsed objects [l^. This model attempts 
to identify such objects from the properties of the linear fields, obtained from the linearization of the 
equations of motion. For our case, this heuristic approach would state that the fraction of matter, 
F{> m), that is enclosed within collapsed objects (here infinitesimally thin shocks, as we consider the 
Burgers equation in the inviscid limit) of mass larger than m, with m = poq, is given by the probability 
that, choosing a Lagrangian point qc at random, the linear-theory Eulerian relative distance xl at time 
t between the particles qc + q/2 and qc — q/2 vanishes. (For a 3-dimensional Universe one considers the 
probability the a sphere of mass m centered on qc has collapsed to a point.) In terms of dimensionless 
variables this reads as 

F^^{> M) = Pq{Xl < 0) at M = Q, with Pq{Xl) = y= , (155) 

where refers to the distribution obtained by linear theory, where particles always keep their initial 
velocity and shocks are discarded. This gives 

FPS(>M)= / dy^. (156) 

As usual, Eq. (|156p implies F^^{> 0) = 1/2, which means that only half of the mass would be within 
collapsed structures. Therefore, it is customary to multiply this by a somewhat ad- hoc factor 2 (ssj . 
Thus, differentiating with respect to M, the standard Press-Schechter recipe gives in our case the mass 
function 

2F^^{>M)^ / dMMN^^{M), whence N^^ (M) ^ ^ Ar^/^ e-^^' . (157) 
Jm 

Therefore, we find that for the 1-D Burgers dynamics with Brownian initial velocity the Press-Schechter 
ansatz happens to give the exact mass function (|154p . The agreement of the Press-Schechter mass 
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function at both small and large masses for the one-dimcnsional case was already noticed in |,51i] . for 
more general power-law initial velocity energy spectra (although there were no exact results available 
at large masses but heuristic predictions). This can be somewhat surprising in view of the many effects 
that could have made the Press-Schechter ansatz fail (especially at the low-mass tail), such as the 
so-called "cloud-in-cloud problem" , associated here with the fact that, even though particles ± g/2, 
evolved according to linear theory, may have not collided yet, on a larger scale ^ > g it may happen 
that particles 1/2 had very large inward velocities and have formed by time t a massive shock 
that includes the smaller scale q. Nevertheless, in the cosmological context, numerical simulations have 
shown that, even though the Press-Schechter mass function is not exact, is usually gives reasonably 
good estimates (e.g., [42|), so that it is still widely used today. It is satisfying to find out that in 
a related dynamical system, it actually happens to coincide with the exact result. This makes the 
reasonable agreement observed in other cases somewhat less surprising than would be expected at first 
sight. As found in [5l|, this also suggests that it could provide a reasonable estimate for the Burgers 
dynamics itself with more general initial conditions. 



8.2 Spatial distribution of shocks 

Finally, from the n— point distributions (|146p we can derive the many-body distributions of shocks, 
far from the origin. Thus, in a fashion similar to Eq. (|152| ). we can relate the trivariate mass function, 
N{Mi,M, M2; X)dMidMdM2dXidX, that counts the probabihty to have a shock of mass Mi within 
[Xi,Xi + dXi], another shock of mass AI2 at distance [X,X + dX], and a mass M in-between both 
shocks, to the three-point conditional shock probability, Pq2i.Q32,Q43j of Eq. (ll47p . This reads as 

/•Q32 roc I'oo /•M1—Q21 

Pq2uQ32,Q.3{X) ^ dM dA/i / dM2N{Mi,M,M2;X) dQi 

"'0 JQ21 JQis Jo 

xe{Qi + Q2i + Q32-Mi~M)e{Mi + M+M2-Qi-Q2i-Q32-Q43), (158) 

where the two Heaviside factors ensure that Q3 and Q4 are within the second shock of mass M2. Then, 
differentiating with respect to Q21 and Q43 yields 

f,^'"^^ ^ / dAf / dM2NiQ2i+Q32 + Q43-M-M2,M,M2;X). (159) 

dQ2ioQi3 Jo Jq43 

Taking the shifted Laplace transform of both quantities, in the form 

iV(Mi, M, M2;X) = /^'" dsidsds2 ^^s,-i)M,+(s-i}M+(s2-i)M2 ^(^^^ x)^ (160) 

and using Eq. (|147p . we obtain 

(V^i + vs)(vs + V«2) 

If we now consider the multiplicity of shocks at positions Xi and X2, independently of the mass 
M in-between, we integrate over M and s, which gives the bivariate mass function of shocks separated 
by a distance X: 

N{M,,M2;X) = «J^(,^.i)M,+(.2-i)M2 4a-^i)(l--^2^ ^ N{M,)N{M2), (162) 

J-^co (271-1)^ (yiT+l)(l + 7iJ) 

where we used Eqs. (|154p and (|14ip to recognize the product N{Mi)N{M2)- Therefore, we find that 
the bivariate mass function ^(toi, m2; x) does not depend on the inter-shock distance x and merely 
factorizes as n{mi) x n{m2) (far from the origin xi =0). 

Thus, shocks are not correlated and there is no bias: knowing that there is a shock of mass mi at 
position xi does not bias in any way the shock multiplicity at position X2 = xi + x. We can note that 
this is consistent with the fact that the densities within separate regions are uncorrelated, as seen in 
Eas. (|114p and (jllSp . In fact, from Eg. piSp . which shows that the density fields over separate regions 
are completely independent, we can see that this must extend to all n— point shock mass functions, 
thus shocks at differeiit positions are uncorrelated. The fact that shocks form a Poisson point process 
was also obtained in [9|| for the case of one-sided Brownian velocity. 
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9 Conclusion 

We have shown in this paper how to derive the equal-time statistical properties of the solution of the 
Burgers equation with Brownian initial velocity, using the transition kernel associated with Brownian 
particles over a parabolic absorbing barrier. This initial velocity field is not homogeneous, as the initial 
velocity is Gaussian with a variance (vq) oc at distance |a;|. However, it has homogeneous increments. 
Then, although the one-point distributions, Px{v) and Px{q), of the velocity v and initial Lagrangian 
position q, depend on the position x, the two-point distributions exactly factorize provided that all 
spatial coordinates remain on the same side of x = 0, such as Pxi,x2i^iT'i'2) = Pxi{vi)Px2ii'^'2i) ^^'^ 
Pxi.x2 {ill 92) = Pxi {qi)Px2i (921) for > 0. A similar factorization holds for higher-order distributions. 
This agrees with the results that were obtained for the one-sided Brownian initial velocity by ^ . In the 
limit where we are far from the origin, this implies that we recover the invariance through translations 
for the distributions of velocity and Lagrangian increments. Then, we have focussed on the properties 
of the system in this limit, where many simple explicit results can be derived. 

As expected, we have found that on large scales, or at early times, all statistical properties converge 
to the Gaussian distributions set by the initial conditions (the nonlinear evolution being subdominant). 
On small scales, or at late times, the distributions of the velocity increment, Lagrangian increment, or 
mean overdensity r/ within a region of length x, increasingly depart from the Gaussian. They exhibit 
widely separated exponential cutoffs, of the form and e~^^^ , while a power law r]~^^^ develops in 
the intermediate range. However, we find that the variance of these distributions remains unchanged by 
the nonlinear dynamics, that is, it is equal to the value that would be obtained by discarding collisions 
and shocks, and letting particles cross each other and always keep their initial velocity. In particular, 
the second-order velocity structure function and its energy spectrum do not evolve with time, while 
the density correlation remains a Dirac function with an amplitude that grows as t^. 

In fact, the densities within non-overlapping regions are uncorrelated and, at any order, the n— point 
connected density correlation can be written as a product of n — 1 Dirac factors: it is non-vanishing 
only when all points coincide. Then, it can also be written as a product of n — 1 two-point correlations, 
that connect the n points, with a constant amplitude that happens to be the number of heap ordered 
trees. This allows a combinatorial interpretation that is similar to the hierarchical tree models that 
were devised as a phenomenological tool in the cosmological context, with the difference that in our 
case we must consider ordered trees (note that in the present 1-D system, where particles do not 
cross, it is meaningful to order particles by their positions so that the concept of ordering appears 
rather natural). Then, the cumulants of the overdensity exactly scale as (77")c esc ^\ with an 

amplitude that is independent of time and scale. Thus, they happen to exactly satisfy the so-called 
"stable-clustering ansatz". In fact, the density cumulant-generating function remains exactly equal to 
the one obtained at tree-order from a perturbative approach (which breaks down beyond leading order 
as next-to- leading corrections actually are divergent, which signals the need for a non-perturbative 
method that takes into account shocks). 

We have also studied the Lagrangian displacement field, associated with a Lagrangian description of 
the dynamics. In the limit where we are far from the origin, we find that it satisfies a similar factorization 
and recovers the invariance through translations for the distributions of relative displacements. This 
also allows us to derive the properties of shocks. Thus, in agreement with previous works, we find that 
all of the mass is enclosed within shocks, at any time t > 0, and that the shock mass function has 
the very simple expression (|154p . It agrees with the asymptotic behaviors that were already obtained 
through analytical means or numerical simulations in or [sH and the exact result of Q . Finally, 

shocks are not correlated as the bivariate multiplicity function n(mi, TO2; x), that counts shocks of mass 
TOi and m2 separated by the distance x, factorizes as n(mi) x 71(7712), in agreement with the same lack 
of correlation obtained for the density field. 

Thus, the equal-time statistical properties of the Burgers dynamics with Brownian initial velocity 
are remarkably simple. It appears that the nonlinear dynamics preserves some properties of the initial 
fields (e.g. the second-order structure functions, the independence and homogeneity of velocity incre- 
ments and of the densities in separate domains) and that simple explicit expressions can be derived in 
the limit where we are far from the origin (or on the right side for one-sided initial conditions). 

At finite distance from the origin, in addition to the quantities given here we need the two-point 
distribution associated with the case where the two particles are on different sides of the origin x — 0. 
Although we can obtain explicit expressions by the method presented in this article, this leads to mul- 
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tidimensional integrals that do not seem to greatly simplify (although we obtained simple expressions 
for the one-point distributions). However, for practical purposes, one is mostly interested in the limit 
where we are far from the origin and homogeneity is recovered. From a physical point of view, the 
initial conditions ^ are meant to represent a system with homogeneous velocity increments, which 
scale as {{Avo)"^) oc \Aq\ as in^ over a finite range, and an energy spectrum Eo{k) oc as in (fTO|) 
over the range of interest. Thus, in practice there would be an infrared cutoff. A, below which i?o(^) 
would grow more slowly than 1/k so that the velocity field is actually homogeneous (in an experimental 
setup there would actually be a finite lower wavenumber, set by the size of the box, and homogeneity 
would only apply far from the boundaries). Then, the initial conditions ([5]) studied in this article can 
be viewed as a convenient mathematical device to represent such a system, with the understanding 
that the special role played by the origin is a mathematical artifact and that only the properties far 
from the origin are meaningful in the physical sense described above. Note that this identification is 
possible because small scales are not strongly coupled to large scales, in agreement with the fact that 
over large scales we recover the initial fields and no strong correlations develop. 

To put this study in a broader context, it may be useful to recall here the main properties of 
"decaying Burgers turbulence" for more general Gaussian initial conditions, ft is customary to study 
the Burgers dynamics |T]) for power-law energy spectra, Eo(k) oc fc" (here we focussed on the case 
n = —2, see Ea. (|TO)) ). Indeed, at late times the asymptotic statistical properties of the velocity field 
no longer depend on the details of the high-fc spectrum (assuming a strong enough falloff ) nor on the 
precise value of the viscosity as a self-similar evolution develops [23l. Issl]. Then, depending on the 
exponent n, the integral scale of turbulence, L{t), which measures the typical distance between shocks 
and the correlation length, and the shock and velocity probability distributions show the following 
behaviors. 

For —3 < n < —1 (which includes the case n = —2 studied in this article, associated with a 
Brownian initial velocity), the initial velocity is not homogeneous but it has homogeneous increments, 
while for — 1 < n < 1 (which includes the case n = associated with a white- noise initial velocity), 
the initial velocity itself is homogeneous. In both cases, the integral scale grows as L(t) ~ ^^/(n+a)^ 
the tails of the cumulative shock distribution and velocity distribution satisfy ln[n(> to)] ~ — to'^^"^, 
ln[ri(> \v\)] ~ — for to oo,\v\ oo, see [U [s^]. At low wavenumbers, below 1/L{t), 
the energy spectrum keeps its initial form, E(k, t) oc fc", whereas at high wavenumbers it shows 
the universal law, E{k,t) oc fc~^, due to shocks [13, [1^. The preservation of the large-scale part, 
E{k, t) oc /c", is associated with the "principle of permanence of large eddies" [13]. Physically, this means 
that, at a given time i, structures of size larger than L{t) have not had time to be strongly distorted 
by the dynamics (in agreement with the simple scaling argument tat,^{x) ^ x which gives x ^ L(t)). 
In particular, not only statistical properties but each random realization is stable against small-scale 
perturbations [1, HI] ■ Then, the tails of the shock and velocity distributions can be understood from 
the initial velocity field. Thus, the velocity difference between the left and right boundaries of a shock 
of mass TO = poq is q/t^ which leads to a probability '-^ e~*^'/'^ /•^^a'^i) ~ e~™ /* (where we did 
not write constants in the exponent) [5l|. We can check that these properties agree with the results 
derived for n — —2 in this paper. 

For 1 < n < 2 the system shows a more complex behavior, since there are three scaling regions for 
the energy spectrum: first a /c" region at very low wavenumbers, below ks{t) ~ ^-1/2(2-")^ next a k'^ 
region between fcs(t) and fci(t) ^ t^^/"^, and finally the standard fc~^ region above khit) [13] ■ Therefore 
the evolution is no longer self-similar. For n > 2 the fc" region disappears (it gives subdominant 
corrections) and the leading-order evolution is again self-similar but independent of n [27, 25]. 

We can hope that the exact results presented in this article for the case of Brownian initial velocity 
could serve as a useful benchmark to test approximation schemes which could be devised to handle 
other initial conditions where no exact results are available. In particular, in the cosmological context, 
the Zcldovich approximation, which corresponds to removing the diffusive term altogether, has alr eady 
been used to test for instance field-theoretic methods that attempt to resum perturbative series [47| . 




The Burgers equation in the inviscid limit, which corresponds to the more efficient adhesion model, 
might also be used for such purposes. In a similar fashion, the general properties of the Burgers 
dynamics (associated with shocks) have already been used to test approximation schemes devised for 




the study of turbulence [13] . 
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Finally, we note that the method described in this article could also be applied to different-time 
statistics, where the parabolas used in the geometrical interpretations would now have different cur- 
vatures. However, we leave such studies for future works. 



A Some properties of the Airy functions 



We recall here some properties of the Airy functions Ai(a;) and Bi(a;) that are used repeatedly in the calculations 
presented in this article. These two Airy functions are two linearly independent solutions to the second-order 
differential equation y"{x) — xy{x) = 0. The first one, Ai(3c), is the only solution that vanishes at both ends, 
X —* ±cx3, whereas Bi(a;) grows to infinity at a; — > -|-cx3 [J. Both are entire functions and they are related 
through 

Bi(a;) = e'"/®Ai(e'2"''3a;) + e-'^/^AiCe-'^'/^), 
while their Wronskian is constant and given by 



Ai(a;)Bi'(x) - Ai'(a;)Bi(a;) = - 



(163) 
(164) 



We also have the integral representation 



(165) 



At a; = we have 



and for Ixl co: 
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(166) 

(167) 
(168) 
(169) 
(170) 



For |Arg(a;)| < 2n/3, the Airy function can also be expressed in terms of the modified Bessel function of the 
second kind Kv as 



|Arg(.)| < f : Ai(x) = iyf^-i/3 {l^'^') , Ai'(x) = -^K,,, (jx'^') . (171) 



Four useful integrals, that may be obtained from the integral representation (|165p . are [Sj 



da;e°^Ai(a;) = e" whence / da; e"^a::Ai(2;) = Q^e 



2 ci3/3 
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(175) 



that could be obtained from the property Ai"(a;) — a;Ai(a;). 
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Finally, using the property [22 



and the relations (I17ip . we obtain 
p > 0, a + /3>0: 



/3f3x 
V 2 



2/3' 



= _L 3-1/6^2/3 



r(i/ + |) \^ 3 6 6 a + p 



(176) 



(177) 



iy>0, a + /?>0 



2/3' 



^ _1 31/6^4/3 - 
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r(i/+|) V 3 6 6 a + 



(178) 



B Half-range expansion and useful integrals 



We show in this appendix how to obtain the solution (|34|) to the half-range expansion problem p2|) -(|33 | ). The 
same method also allows us to derive other useful identities that we need to perform the calculations presented 
in this article. Thus, we consider the function f{p) of the complex variable p defined by 



s > 0, u > 0, lArg(p)l < TT : f{p) 



Ai 



1/3 I 
p u + ■ 



,2/3 



(179) 



This function is regular over th e com plex plane except for a branch cut along the negative real axis. Moreover, 
from the asymptotic behavior p67p of the Airy function we obtain 



u>0: /(p) 



-1/4 



2V^ 



-1/4 



2„3/2„l/2 



e 3' 



as \p\ — > 00 with |Arg(p)| < n, 



and 



u = 0: /(p)~Ai(0)p ^^'^ as \p\ ^ 00 with lArg(p)l < vr. 
Next, we introduce the general integral Fk.i{vi, .., Vk', Ai, .., A^) defined by 

f c+ioo 

c > max{ Pi } , 



dp 



fip) 



with 



(180) 
(181) 

(182) 
(183) 

If = or £ — one of the products in Ea. p82|) is removed and replaced by a factor 1. If fc = the contour 

in Ea. (|182|l again runs to the rig ht of all si ngularities, that is c > 0. 

Then, from the asymptotics p8Up - p81|l . we can see that we can push the contour in Ea. (ll82|l to the right, 
as c +00, which shows that Fk,e = (using k + £ > 1). On the other hand, by pushing the contour to the 
left, using again the asymptotics (|180p - (|18ip . we can see that Fk,e is the sum of the k residues at p = vf and 
of the contribution associated with the branch cut along the negative real axis. This yields 



with the conditions 



Ui > 0, / for i / i'; Aj > 0, Aj / Xf for j / j'; k + £ > 1. 
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[1,(P + A?) 



(184) 



where C is the anticlockwise Hankel contour that bends around the negative real axis. Then, pushing the 
contour towards both sides of the negative real axis, making the change of variable p — 
and e ^ 0"^, and using [J| 



Ai [e±'^^/^x] 



1 



±i7r/3 



[Ai(a;)TiBi(i)], 



'[i zb ie with /i > 
(185) 
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as well as the Sokhatsky-Weierstrass theorem, written here in concise form as 

1 



lim 



-.0+ M 



— = p.v.— !— - i^sip - v), 



we obtain 



(186) 
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(187) 



Here, the symbol p.v. stands for the Cauchy principal value and must be understood with respect to fi'^ (rather 
than jj,). That is, the integrals are regularized by cutting around each pole Xj the interval [(U_,/x+], with 
fj,± = Aj ± e, which is symmetric in terms of /i^ around A^, and taking the limit e — > 0"''. If f = the integral 
is regular and there is no need to introduce the principal value. In particular, the case k = 1 and £ = yields, 
with ;/ > 0, 



u> 0: 



d^ jM^g- 2.3/2, 
27r fj? + 1/3 
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-fiu H ^ 
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(188) 



This implies that (j}s,u{r, u) defined by Eq. (|34|) is the solution of the form (|32|l that satisfies the constraint (|33p . 
Note that the restriction to it > is essential. For instance, as Arg(it) grows from 0, the function f(p) displays 
an exponential growth for tt — 3Arg(u) < Arg(p) < tt and we can no longer bend the integration contour onto 
the negative real axis. 

Alternatively, as in [Til ], we can obtain Ea. (ll88p from the analysis of 31] and ^32], or of 28j, who studied 
several problems associated with the Klein-Kramers equation, by taking the limit of zero friction. However, 
these problems lead to discrete spectra and require a sophisticated analysis that involves infinite products to 
handle the poles associated with all eigenvalues. 

Finally, making the change k ^ k + 1 and taking the limit Vk+\_ —> 0^ in Ea. (|187|) gives the useful identity 



u> : p.v. 
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'XiU + ■ 



Equation (|189|l now applies to any fc > 0, ^ > 0. In particular, the case k = £ = yields 



(189) 
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(190) 



Note that s can be absorbed in Ea. (|190|l through the change of variables i; = -^su and v = By letting 

m parameters Vi going to zero in a sequential manner, we could derive a series of similar identities for integrals 
of the form of (|189p wit h a p r efact or ^3/2-3m m > 0. However, in this article we do not need to go 

beyond m = 1 as in Eas. (fT89)) -(fT90 )> . 

In a similar fashion, we now consider the function f{p) of the complex variable p defined for Si > 0, Ui > 0, 

by 
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(191) 



It has a branch cut along the negative real axis and it grows as most as p^^^^ for \p\ ~* oo and < Arg(p) < n. 
Then, in a manner similar to Eq. (|182p . we introduce the integral Fk.tiyi, .., I'k', Ai, .., \i) defined by 



Fk,i(i'i', Aj) 



2- ntib-'^f) n=i(p+A= 



with c > 0, 



(192) 



with the conditions 



Vi > 0, Ui ^ Vii for i / i'\ Xj > 0, Aj / Xy for j ^ j'\ k + £> 2. 



(193) 
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Note that the integration contour is now parallel to the real axis, in the upper half-plane, and that the 
asymptotic behavior of /(p) for large \p\, with < Arg(p) < vr, now requires k + £ > 2. Pushing the contour 
upward, as c — > +00, we can see that Fk,t = 0. Next, pushing the contour towards the real axis, by making the 
change of variable p = + ie with /i > and e ^ 0^, and using Eq. (|185|) . that also yields 



Ai' [e^'^^^^x] = ieTW3 [Ai'(a;)TiBi'(x)] 



(194) 



as well as Ea. (|186p . we obtain after taking the real part 
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(195) 



Taking the imaginary part gives another identity, that involves t he in tegral over /i of products Bi 'Ai and Ai'Bi, 
which we do not need for the present calculations. Again, in Eq. (|195|) the Cauchy principal value is understood 
with respect to fj,^ . Exchanging the derivative in expression (|191|l . we obtain an identity similar to Ea. (ll95|l 
where the derivatives are exchanged; 
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i=l 

Next, making the change k ^ k + 1 and taking the limit Vk+i — > gives, for k + £ > 1, 
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(197) 



Here we combined both Equations (|195|l - p96p . and in each bracket the second product, such as AiAi', is equal 
to the first product where we exchange the derivative. Next, for the case ui — U2 = 0, we again make the 
change k ^ k + 1 and take the limit Vk+i 00. This yields for any fc > 0, ^ > 0, 
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where we used the Wronskian property (|164|l and 5k+e,o is the Kronecker symboL In particular, for k = £ = 
we obtain 



Ai'(^)Ai(^)-AiAi' 



(199) 



Since the integral is convergent ai fi — there is no need to use the principal value. 
C Computation of the two-point distribution 

We present here the computation of the two-point distribution (|82l) from the two contributions and 

described in Figs. |4]and[5l 

Let us first consider the contribution p^. In a fashion similar to Ea. (|55p . using the Markovian character of 
the process q t-^ {tp, v}, it reads as 

Pxi,x2{0 < q'l < qi,ci;q'2> q2)dci = lim di>-dv-Kxi,ci{0,0,0;q-,ip-,V-) 

i3±— >±oo / 



X J dV'idui[Jfa;i,ci(0,0,0;gi,?/>i,i;i) - if:ci,ci+dci(0,0,0;(ji,i/>i,-ui)] 

X |y dil!+dv+Kxi,c-i{qi,iJi,vi;q+,tlj+,v+) — J dili2dv2Kxi,ci{qi,^pi,Vi;q2,i^2,V2) 

X d^l;+dv+Kx2,c2{q2,i'2,V2;q+,ij+,v+)}. (200) 



Then, we recognize p^-^ {0 < q'l < qi) in the contribution associated with the first term in the bracket, and we 
can write 

Pxi,x2{0 <q'i< qi;q'2 > q2) =p^i(o < gl < gi) -pJ^lxsIO < gi < gi;g2 > 52), (201) 

where we introduced the remaining part 

^11.12 (0 - 91 — 9i' "^1' 92 > g2)dci = lim / d^_dii_d'(/)idwid-(/'2d«2d'(/)+du+ 
X ,ci (0, 0, 0; q-,tp-,v-) [K^i ,ci (0, 0, 0; gi , V"! , wi ) - -ftT^i ,ci+dci (0, 0, 0; gi , i/^i , ui )] 

xKxi,cAqiy'4^i,vi;q2,i'2,v2)Kx2,c2{q2,Tp2,v2;q+,'ip+,v+). (202) 

In Ea. (|20ip we have also integrated p^ over ci. 

We can note that Pxi,x2 ^^d Pxi,x2 satisfy the following boundary conditions. First, taking the derivative 
with respect to gi to obtain the probability density Pxi.x2{qij 92 ^ 92), we have 



lim^p^j,^2(9i;92 > 92) ^ Pxi(gi), whence lim^ p> (gi; ga > g2) = 0. (203) 



Indeed, following the discussion at the beginning of section [STTI if ^o(g) is tangent to 1^x1,01 at gi, the second 
parabola 'Px2,c2 can only cross the first one at a point g, > gi (otherwise, t/jq being continuous, if we had 
g, = gi the curve V'o(g) would go below 1^x2,02 just beyond gi). Then, all curves tangent to Vxi,cx at gi satisfy 
both properties g* > gi and g2 > gi, so that they are all included in the contribution P^-^^X2 ^ take the 
limit g2 q^ and we must recover Pxi{qi), as stated in (|203p . Second, for large g2 we obviously have the 
asymptotics 

lim P^i,x2(0 < gi < gi;g2 > g2) = 0, lim p^i,^:^ (0 < gl < gi; g2 > 172) = Pr^i (gi). (204) 

q2— * + q2— * + oo 

This latter constraint can be d irec tly checked on Eq. (|202p . 
Using the transformations (|17|l and (|20p . we obtain 

^21 _ <?2 f 

Pxi,x2{0 < q'l < qi;q'2> q2) ^ e ^ driduidr2du2dr3Haoir3,ui) 

xZ\(gi; rs, — ui; n, ui)G(g2 — gi; ri, ui; r2, U2 + U2i)Hao{r2,U2), (205) 
where we introduced as in (|57p the quantities 

^'^V"nT' ^21 =62 -Ml- (206) 
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Nex t, taking the derivative with respect to qi, using the backward equation p7|l for G and the forward equation 
(|22|) that is also satisfied by A, we obtain a total differential over ri, which only leaves the boundary term at 
n = 0: 

, "21 12 f 

Pxi.x2ili'^<l2 > <12) = e '' ^ / dMidr2du2dr3 //■oo(r3,ui)uiZi((7i; ra, -Mi; 0,-ui) 

xG{q2 — qi;0,ui;r2,U2 + U2i)Hcx,{r2,U2). (207) 

We can check from the explicit expressions of A and G that the integrations by parts leading to Eq. (|207p are 
valid. The expression (|207|) clearly s atis fies the property (|203[) . Indeed, for q2 — > q^ the factor G implies r2 — > 0, 
using the firs t bo undary condition (|23|) . which in turns leads to U2 + U21 < because of the secon d b oundary 
condition in (|23p . However, the last factor Hoc also implies U2 > 0, using the boundary conditi on (139 p applied 
to Hoc- Since U21 > b oth constraints on U2 cannot be simultaneously satisfied which leads to p03|) . This can 
also be checked on Eq. (|2U7p using the explicit expressions of A and Hao ■ 

Then, using the ex plicit expressions of G, A, and Hoo, and the results of Appendices |X] and [B] it is possible 
to greatly simplify Ea. (|207p . Indeed, the integrals over are immediate (they only involve factors of the form 

e"" whereas the integral over ui can be transformed using Eq. (|175p . Next, integrals over Ui are typically 
split over Ui < and Ui > 0, and each factor of the form Ai(— / iMj + s/ y?) with Ui > 0, or Ai{fj,Ui + s/fj,^) with 
Ui < 0, can be integrated over fj, using the results p87p - p90p of Appendix [B] This leads to products of the 
form Ai{iiUi + si//i^)Ai(/xUi + S2//i^) that can be integrated over Ui using the primitive (|63p . that also extends 
to the s econ d Airy function Bi. Then, these terms can be further simplified using the Wronskian (|164p and the 
results p95p - (|199p of Appendix[Bl We eventually obtain in terms of dimensionless variables 

/ + joo , , 
HgH^g(si-i)Qi+{.2-i)Q2ij(,^^2Xi)/(s2,2X2i), (208) 



with 



>i>0, O2i=O2-Qi>0, and /(s,2X) = — !— e-(§^'^'-^+3)2-^/(''-i). (209) 

s — 1 



We can check that Ea. (|208p agrees with both constraints (|203p - (|204p . Then, taking the derivative with respect 
to Q2 we obtain the full probability density associated with g, > q2 as 

dgd^ e(^-^)'3^+(^-^'«- J(.i, 2X0(.2 - 1)I(.2, 2X21), (210) 

-zoo ^ ^ 



since the derivative with respect to Q2 of the first term in the right hand side of (|20ip vanishes. 

We now consider the second contribution, , associated with the intersecti on g , between both parabolas, 
1^x1,01 and 7^X2, C2, being in the range gi < g, < 52- Proceeding as for (f55)) and pOOp it reads as 

Pxi.x2{0 <<li< gi,ci;g2 > g2,C2)dcidc2 = lim / dip-dv-Kxi,ci{0,0,0;q-,ii-,V-) 



X j d7/)idwi[_K'a;i,ci(0,0,0;gi, V)!,"!) - -ft'i^i,ci+dci(0,0,0;gi,?/)i,iii)] 

X j Ail)tdVtKxi,ci{q\,ipi,vi\qt,il)t,Vt) j Aii2dv2Kx2,c2{q*,i'*,v*\q2,ip2,V2) 

X / '^'4'+'^V+[Kx2,C2{q2,'4>2,V2;q+,'4!+,V+) - K^2^^2+ic2{<l2:i>2,V2;q+..'4>+,V+)], (211) 



which we must integrate over both ci and C2. We can note that it satisfies the boundary conditions 

lim P^i,^2('3'i;'?2 > ga) 0, and lim p^,,^2 (0 ^ ^ 92 > ^2) = 0. (212) 

52->g+ 92^+00 

Using again Eas. (|17p . (|20p . we obtain 

"21 _ 92 r 

Pxi,x2{0 < q'l < <li;q'2 > 12) = e '' I dridMidr,du*dr2du2dr3dr4-H'oo(r3,Mi) 

dH 

xZ\(gi; rs, — Mi; n, ui)G(g. — gi; ri. Mi; r*, Mt)G(g2 — g.;r*,Mt — M2i;r-2,M2) „ °° (r2,M2). (213) 

ar2 



40 



Next, taking again the derivative with respect to qi and using the forward and backward equations (f22|l . p7|l . 
gives 

"21 12 r 

Pxj_,x2i<li;<l2 > q2) = e -I ^ duidr,du,dr2du2dr3dr4Hao{r3,ui)uiA{q-i;r3,-Ui;0,ui) 

OH 

xG{q, - gi; 0, ui; r*, it*)G((72 - q,;r.^,,Ut - U2i;r2,u2) {r2,U2). (214) 



As for the derivation of Ea. (|208|) . usin g the exphcit expressions of G, A, and Hoo, and the results of Appendices 
lAlandlBl as well as the property (I28|l . we obtain 

ClSidSClS2 siQi + s(Q,-Qi)+S2(Q2-Q.) 



xJ(si, 2X1)— !—-[L(l,s; 2X21) -L(s2,s; 2X21)], (215) 

S2 — 1 

with 

L(si,S2;2X) = ^^^'^'^^ i(s,s;2X)= 2^e-^^^. (216) 

We now have three inverse Laplace trans forms because of t he three terms AGG in Eq. (|214|) . T he in tegration 
over Q, is associated with in Ea. (|214p and C2 in Ea. (|211|) (g* being related to C2 through Eq.(|8T}). It gives 
a factor [e^''^^ — e"^''^^'^]/ [s — S2). Then, choosing for instance a contour such that K(s2) > 3f?(s) > 1, we can 
integrate the first term over S2, which gives zero by pushing the contour to the right, 5R(s2) +00, and the 
second term over s, which gives the contribution associated with the pole at s = 82- This yields 



Pi„X2(Qi,Q2 > Q2) = 2X21 e^^- 



g(si-l)Ql+{s2-l)Q21 



X J(si, 2X1)— !— -[L(l, S2; 2X21) - L{s2, S2; 2X21)], (217) 

S2 — 1 



We can check that Ea. (|217p agrees with the constraints (|212l) . Then, taking the derivative with respect to Q2 
we obtain the probability density 



dsidS: 



i g(31-l)Ql+(s2-l)Q21 



(27ri)2 

X J(si, 2Xi)[L(s2, S2; 2X21) - L{1, 32-, 2X21)]. (218) 

Finally, combining Eqs. (|210|l and (|218|) we find that two terms cancel out and we are left with the total 
probability density 
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